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1.   INTRODUCTION 

The  dipole  antenna,  one  of  the  most  useful  and  simplest  of  all  an- 
tennas, has  been  and  is  still  an  object  of  interest  to  antenna  theore- 
ticians.  In  all  but  the  elementary  geometrical  shapes,  this  antenna 
poses  a  cryptic  problem.   Recently  the  fat  stub  antenna,  a  dipole  or 
monopole  whose  transverse  dimensions  are  comparable  with  the  wavelength, 
has  aroused  some  interest  as  a  possible  instrument  for  radiating  energy 
over  large  azimuthal  angles.   By  varying  the  shape  and  size  of  the  mono- 
pole  and  by  changing  the  feed  geometry,  wide  variations  in  radiation 
patterns  can  be  obtained.   Electrically  rotated  search  beams  are  another 
possibility  as  will  be  mentioned  later. 

Spherical  antennas  have  received  much  attention  in  the  literature. 
Stratton  and  Chu^-'-''  and  Schelkunof f ^'  have  treated  the  uniformly  excited 
or  rotationally  symmetric  case,  finding  solutions  for  the  pattern  and 
radiation  conductance.   A  spherical  surface  with  external  wires  has  been 
studied  by  Feld'"^'.   DuHamel'^'  has  treated  a  system  of  currents  flowing 
on  the  sphere,  from  the  point  of  view  of  synthesizing  patterns.   A  novel 
continuation  of  Stratton  and  Chu's  work  has  recently  been  made  by  Karr'^)] 
who  considers  a  sphere  split  along  a  circle  of  latitude  not  the  equator, 
and  uniformly  excited.   All  of  these,  however,  yield  far  field  distri- 
butions that  are  invariable  about  the  azimuthal  angle  and  involve  only  a 
simple  set  of  modes. 


(1)   See  numbered  reference  in  the  Bibliography.   All  other  such  numbers 
used  throughout  this  manuscript  also  refer  to  numbered  references  in 
the  Bibliography. 
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It  is  the  purpose  of  this  paper  to  examine  in  detail  the  most  prosaic 
prototype  of  the  fat  stub  antenna,  the  spherical  dipole  with  asymmetric 
feed.   One  may  think  of  this  antenna  as  being  made  by  taking  a  metallic 
sphere  (hollow),  removing  a  narrow  strip  of  metal  around  the  equator  to 
divide  the  sphere  into  two  hemispheres,  and  then  applying  an  electric 
field  in  the  narrow  gap.   If  the  electric  field  is  a  function  of  the 
azimuthal  angle,  then  the  excitation  is  said  to  be  rotationally  unsym- 
metric  or  just  asymmetric. 


2.       THE   GENERAL   ANTENNA   PROBLEM 

2.1      Introduction 

In  attempting  exact   solutions   of  antenna  problems,    Maxwell's   equa- 
tions  in   differential    or   integral    form  are   usually  used   as   a   starting 
point.      The    former   are    for   a   homogeneous   and  isotropic   medium 

m  9E 

curl   E  =  -u.  — ,      curl  H  ■   e  r-.  (2-1) 

The  divergence  equations  have  been  omitted.  These  equations  are  equiva- 
lent to  the  vector  wave  equation  in  E  (or  H)  obtained  by  the  elimination 
of  H  (or  E)  in  Eq.  2-1 

32E 

-  curl  curl  E  +  ^e  — —  =   0  .  (2-2) 

at 

It  is  often  easier  to  deal  with  the  vector  Helmholtz  equation  obtained 
from  Eq.  2=2  by  assuming  a  monochromatic  time  dependence  eJwt: 

curl  curl  E  +  k*E  =  0  (2-3) 

where  k  =  w  [ie    . 

A  vector  in  space  is  determined  by  its  three  scalar  components  and 
clearly  the  problem  would  be  immeasurably  simplified  if  the  solutions  to  the 
scalar  Helmholtz  equation,  about  which  a  great  deal  is  known  could  some- 
how be  used.   Since  Maxwell's  curl  equations  include  div  E  -  0,  E  can  be 
written  as  a  curl,  i.e.,  as  a  solenoidal  (transverse)  vector.   The  latter 
is  determined  by  two  scalars  only  as  the  div  E  =  0  condition  imposes  an 
additional  linear  relationship.   In  certain  cases,  the  vector  Helmholtz 


4 
solutions  are  found  by  taking  the  curl  and  curl  curl,  respectively,  of  a 
suitable  unit  vector  times  the  scalar  Helmholtz  solutions.   The  cases  in 
which  this  is  possible  are  six:   the  six  coordinate  systems  which  are 
listed  below^o).   jn  three  other  cases,  (spheroidal  coordinates  and 
rotational  parabolic  coordinates)  if  the  cp  variation  is  absent,  (m  =  0) 
the  vector  equation  can  be  reduced  to  a  single  scalar  equation  which  can 
be  solved.   For  other  systems,  solutions  to  the  equation  might  be  found 
but  fitting  the  boundary  conditions  is  almost  impossible.   The  separable 
systems  are:   cartesian,  with  any  axis  the  preferred  (longitudinal) 
axis,  spherical  and  conical,  with  the  radius  the  preferred  axis;  circular, 
elliptic,  and  parabolic  cylinder,  with  the  z  axis  preferred.*  Solution  of 
the  vector  equation  is  then  reduced  to  finding  scalar  solutions  and 
operating  on  these  scalar  solutions  to  give  the  vector  solutions.   It  is 
well  known  that  the  scalar  Helmholtz  equation  can  be  solved  in  any  co- 
ordinate system  which  allows  separation  of  the  partial  differential 
equation  into  three  ordinary  differential  equations,  each  depending  on 
only  one  coordinate.   There  are  eleven  such  systems^''  and  the  six  systems 
mentioned  above  are  among  these  eleven. 
2  2  The  Spherical  Antenna 

The  system  of  interest  is  the  spherical  coordinate  system,  since  the 
boundary  conditions  are  to  be  applied  over  a  spherical  surface.   The 
orientation  is  shown  in  Fig.  1. 

If  the  scalar  Helmholtz  equation  is  written  in  the  form 

V2f  ♦  k2f  -  0  or  V2g  +  k2g  -  0,  (2-4) 


Figure  I.  Spherical  Coordinate  System 


two  independent  vector  solutions  of  Eq.  2-3  are  given  by 


(8) 


Mr  =   curl   a   rf,  N_  =   f-  curl   curl   a   rg 

_1  _x  _g    b  _I 


(2-5) 


The  components  of  these  vector  solutions  along  the  coordinate  axes  are: 


N. 


Mr 

"  0 

n(n    + 

D< 

k 

I 

W 


J 91 


N< 


r  sin  9  3qp 


=  JL 


% 


kr  3r39 


Nn 


-  1  9£ 
r  80 


JL 


(2-6) 


kr  sin  0  3r3cp 


For  one  set  of  fields  let  E  =  Mf;  then  H  =  j/r\  Nf  with  n,  a  wu./k  =  -T\jje 
The  other  set  is  taken  as  E  -  N„  and  H  =  j/r\  M„.  These  are  seen  to  be 
transverse  electric  and  transverse  magnetic,  respectively.  The  scalar 
solution  has  two  separation  constants  or  eigenvalues  n  and  m  and  is  of 


the    form 


g(r,9f(p)    =   h(2)(kr)    t^U^P^Ccos   9)    cos 


mcp 


(2-7) 


+  BnmUnmPnm(cos    9>    sin  m(P^    or 
f(r,e,cp)    -   h(,2)(kr)     [CninUninPn,n(coS    9)    cos  mcp 

!"  DnmUnmpnm(cos  9)  sin  m^  > 

where  h^^^kr)  is  a  spherical  Hankel  function  of  the  second  kind,  repre- 
senting an  outward  travelling  wave;  P  m(cos  9)  is  the  associated  Legendre 
polynomial  of  degree  n  and  order  m.   General  references  on  spherical 
Bessel  functions  are  Stratton'  "' ,  and  Erdelyi'  ^' .   For  associated 
Legendre  functions,  see  Hobson^*',  Erdelyi^"^,  MacRobert^^',  and 
Whittaker  and  Watson'^"*)    Excellent  tables  are  given  by  the  National 
Bureau  of  Standards^^,  *•"'..   The  spherical  Hankel  function  is  related  to 
the  more  common  cylindrical  Hankel  function  by 


hn(2)(kr)  .  -2-^Vy2   (kr) 


Ttk] 


(2  8) 


n  n.kr 


That   these    functions    are   simple  may   be   seen    from   the   expressions    for   the 
first    two,    which   are 


h0<2)(x)    =  ?ie-Jx         h,(2)(x)    ■   z2e  jx(1  .  J} 
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Similarly   the   associated  Legendre    function    is    related   to   the  ordinary 
I • "endre    function    by 


Pnm(co.  9)   -  sin"  9  ^  Pn(cOS  6) 

d   (cos   9)m 


(2-9) 


Unra  is  a  normalizing  factor  which  will  be  defined  later. 

The  transverse  components  of  electric  and  magnetic  field,  called  ET, 
ftp,  can  now  be  written  in  terms  of  vector  normal  mode  functions  eK(0,qp) 
and  mode  voltages  Vx(kr): 


h  -  Y_v* 


(kr)  eK(0,<p)  (2-10) 


ftp  -  2_  \   V<kr)  »r  x  eK(0,cp)  ,  (2-11) 


where  K  means  a  summation  over  all  possible  modes.   The  prime  on  V)./(kr) 
indicates  differentiation  with  respect  to  r;  j/^k  V^'Ckr)  can  be  called  a 
mode  current.   The  mode  functions  and  voltages  are  listed  here  for  the 
various  modes. 

TE  even: 

Vn[n(kr)    =  C^rh^^kr),      e(6,<p)    ■  -ar  x  V[UnmPnm(cos   0)    cos  mcp]      (2-12) 

TE  odd: 

Vnm(kr)    =  E^rhn^^kr),      e(9,<p)   ■  -ar  x  V[UnmPnm(cos   0)    sin  m(p]      (2-13) 

TM  even: 

« 

Vnm(kr)    =  ^[rhj',      e(0,(p)   ■  V[UnmPnm(cos   0)   cos  mqp]  (2-14) 

TM  odd: 

Vnm(kr)   =  ^   [rhn]',      e(0,cp)   ■  V[UnmPnm(cos   e>   sin  **]  (2-15) 

where 

frhj'   -  f  [rh^Ckr)]  (2-16) 

or 


It  can  be  shown  that  for  TE  modes, 

V'(kr)  =  *'V '  V(kr) 


(rh^ 


while  for  TM  modes, 


V'(kr) 


-k2rh 


?  V(kr)  6 


2.3  Evaluation  of  the  Mode  Voltages 

From  the  Uniqueness  Theorem^1' '    for  electromagnetic  fields,  the 
exterior  field  is  completely  specified  by  tangential  E  over  the  surface 
of  the  sphere,   Over  the  conducting  portion,  the  tangential  electric 
field  is  zero.   Thus,  specifying  tangential  E  over  the  gap  or  aperture 
uniquely  determines  the  field  outside  the  sphere. 

The  basic  method  here  is  to  assume  a  suitable  distribution  of  Etan 
in  the  aperture,  to  express  the  external  and  internal  fields  in  terms  of 
the  assumed  distribution,  and  then  to  make  the  two  forms  for  Htan  agree 
at  the  gap.   For  a  narrow  gap,  the  actual  distribution  of  electric  field 
(except  for  the  distribution  with  respect  to  qp)  in  the  gap  does  not 
appear  to  affect  the  performance  of  a  dipole  antenna  much.   It  will  be 
assumed  that  the  electric  field  in  the  gap  does  not  vary  with  z,  and 
that  it  has  only  a  z  component.   For  narrow  equatorial  gaps,  Ez  -  Eq. 

Broadly  speaking,  the  external  and  internal  fields  are  found  in 
terms  of  the  gap  field  by  equating  the  assumed  distribution  to  the 
general  series  form  (Eq  2-10),  multiplying  both  sides  by  one  mode, 
integrating,  and  using  the  orthogonality  of  modes  to  reduce.   This 


gives  a  result  of  the  form 

Vq(kr)  -  P^|  (2-17) 

where  q  is  one  mode  of  x.   This  procedure  will  be  followed  in  detail  in 
the  next  few  pages.   Call  the  assumed  gap  field  Eg  =  Eq(cp).   A  convenient 
method  uses  the  Lorentz  Reciprocity  Theorem^  "''.   This  theorem  relates 
two  electromagnetic  fields  E1 ,  Hi  and  E2 ,  H2  of  the  same  frequency  in 
space  by  an  integral  over  the  surface  enclosing  the  source  free  region,  G. 

/GEi  *  H2-n  dS  =  ^G?2  x  Ht'n  dS  .  (2-18) 

Here  the  region  is  taken  as  that  enclosed  by  the  antenna  surface  of 
radius  a   and  a  sphere  of  large  radius  r,  r-*00.   The  fields  must  satisfy 
the  radiation  condition^^'  at  infinity 

Lim  r[n  *  H  +  1  E]  =  0  (2-19) 

r-oo  -  -  n  - 

Let  the  antenna  surface  be  Sx  and  the  surface  of  large  radius  be  S2. 
Then  the  integral  of  (Ei  *  H  -  E2  x  Hi )  over  Si_  must  equal  that  over  S2. 
Now  let  r  vary.   Since  the  integral  over  Si  has  not  changed,  that  over  S2 
must  also  be  a  constant.   If  the  limit  of  the  integral  over  S2  is  taken  as 
r-°°,  and  the  radiation  condition  used,  zero  results.   Thus,  the  constant 
value  is  zero.   This  leaves  the  integral  over  the  surface  of  the  antenna 
plus  gap;  however,  due  to  the  perfect  conductivity  of  the  sphere  the 
integral  over  the  gap  surface,  S,  alone  remains. 

Let  the  field  corresponding  to  the  assumed  distribution  of  tangential 
E  be  subscripted  1,  whereas  that  subscripted  2  will  be  one  TE  (or  TM) 
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mode.   The  single  mode  field  may  be  further  split  into  even  and  odd  terms 
because  the  sin  mqp,  cos  mqp  factors  are  orthogonal.   H±  is  unknown:  thus, 
the  general  series  expression  (2-11)  must  be  used^  u'.   Four  cases  are  to 
be  considered  and  each  case  should  yield  one  set  of  mode  voltages.   One 
of  these  cases  will  be  examined  in  detail. 
Case    1  E2 ,  H2  are  TM  even. 

h   =   Vmn£nm>       H2  =  i  V^  ar  x  enm  (2-20) 

with  the  fields  being  evaluated  at  r  =  a.   The  left  hand  side  of  Eq.  2-18 
is 

2 

which  reduces  to 

.   2 

^   "   ^  Vnm  CJe*pQ  ^'fnm   sin   ed9d^  (2-21) 

The    right   hand   side   of  Eq      2    18    is 

R"5  =   a2C!o    Vnm4  !nm  x  Z    Z    (ar  X  !ji^r  Vji   sin   edGdcp 


which  simplifies  to 

RHS    "    —  Vnm  CC   IIVji   £jx    £nm  sin   9d0d*  .  (2-22) 

nk         j   j 

It  is  now  possible  to  interchange  the  order  of  summation  and  inte- 
gration since  the  series  of  Eq.  2-22  is  uniformly  convergent  in  9  and  qp. 


11 

It  will   be  shown   that   the  vector  mode    functions   are  orthogonal: 
fo    ;o^i°fnmdS  =   r  4    ;oV[Un  Pnmcos  "»]  "VtUj^cos   i<p]sin   9d0d<p 

=   42"-CCUninpnin'UjiPji'cos  mcp  cos    icp   +  -M-  Unn,PnmU, *p, isin  mcp  sin  i<p] 

sin  6 

sin   9d9dcp 
■   6nj6mi  2*n(n  +  1)  ,  (2_23) 


em 


This  can  be  seen  from  the  expressions  for  orthogonality  of  the  trigono- 
metric and  associated  Legendre  functions '°' ,  which 


are 


r2K  «  •  2K 

Jq    cos  mqp  cos   lcpdqp   =   *±  6^,    J      sin  m<p  sin  itpdqp  =   «(l-6m0)fim1  (2-24) 

m 

where  .         [1      if  i    =   i 

0     if  i   f  m        is    the  Kronecker  delta'^JJ^ 


m 


em 


2     m  =   1,2,3, 

1     m  =   0  is   the  Neumann  number^^',    and 


Tt   dP  m  dP  m  2 

!° ^  Sr +  rize  W]'»  9de  ■  6-J  ^fr^r^     <2-"> 


The  normalizing  factor  U_m  is  defined  as  the  square  root  of  the  normal 
ization  factor  for  the  associated  Legendre  polynomials  and  is 


Unm 


(2n  +  l)(n  -m)!  (2_26) 

2  (n  +  m) ! 


12 
This  factor  has  been  used  because  of  convenience  in  the  computer  sub- 
routine for  calculating  the  polynomials;  this  code  will  be  discussed 
later.   The  result  of  these  manipulations  is 

RHS  =  —  V  V  '  2nn(n  +  1)  (2-27) 

-nh      nnrnm     P  K£-  *•*■  ' 

MK  fcm 

with   the    final    form   for   the  mode   voltages 

Vnm(ka)    =   2rtn(nmr  1)   ^W^  !e'!nmsin  0d^   •  (2-28) 


Case    2  The  remaining  cases  are  identical  in  result  (Eq.  2-28).   The 

differences  appear  when  the  values  for  the  mode  functions  and  voltages 
from  Eqs=  2-12  through  2-15  are  inserted  in  Eq.  2-28. 

To  recapitulate,  the  electric  and  magnetic  fields  in  the  exterior 
region  were  written  as  a  summation  of  mode  functions  and  voltages  in 
Eqs .  2-10,11,  and  the  mode  voltages  have  been  found  in  terms  of  an 
assumed  distribution  of  tangential  E  in  the  gap,  Eq.  2-28.   Next,  the 
interior  fields  will  be  handled  in  a  similar  manner. 


3.   THE  RADIAL  LINE 

3  I  The  Radial  Line  Exciter 

The  spherical  antenna  of  the  last  chapter  was  excited  by  an  electric 
field  distribution  in  the  equatorial  gap.   It  is  the  purpose  of  this  sec- 
tion to  consider  methods  of  producing  various  voltage  distributions  in  the 
gap.   Since  it  is  usually  desirable  to  feed  practical  antennas  with  a 
coaxial  cable  whose  diameter  is  small  with  respect  to  wavelength,  some 
sort  of  transmission  line  or  waveguide  is  needed  to  transform  from  the 
cable  to  the  large  diameter  gap.   A  number  of  cables  might  be  used  to  feed 
the  gap  directly  but  this  is  an  unwieldy  solution.   Eiconical  and  parallel 
plate  transmission  lines  appear  to  be  two  feasible  types.   The  biconical 
line  is  simple  of  solution  only  as  long  as  the  apex  is  central  with  the 
sphere.   For  this  reason  the  parallel  plate  line  was  chosen  to  handle 
more  complicated  feed  geometries.   The  parallel  plate  line  consists  of  two 
perfectly  conducting  parallel  plates,  spaced  a  distance  apart  equal  to  the 
gap  distance,  located  interior  to  the  sphere  and  connected  to  the  sphere 
at  the  gap.   Thus  any  radial  component  of  current  flowing  on  the  parallel 
plate  line  can  flow  through  the  gap  out  onto  the  surface  of  the  split 
sphere.   If  the  line  is  excited  at  one  point,  it  may  be  called  a  radial 
transmission  line.   A  cross  section  of  sphere  and  radial  line  is  depicted 
in  Fig.  2. 
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Figure  2=  Cross  Section  of  Antenna  with  Radial  Line 

3.2  The  Field  Solutions 

Consider  a  radial  transmission  line  with  the  orientation  of  Fig.  3. 
Using  cylindrical  coordinates  and  the  techniques  of  the  previous  chapter, 
a  vector  wave  equation  solution  is  constructed  from  the  scalar  wave  equa 
tion  solution.   Let  the  scalar  solutions  be  f  and  g   The  most  convenient 
vector  solutions  are  then  given  by  the  following: 


*f 


curl  a  f 
-z 


N  -   1/k  curl  curl  a  g 

-g  -■ 


(3  1) 


Mr   -    1/k  curl  N, 
-f  -f 


The  components  of  these  are 


1/  9f 
1/r  — 

dtp 


9f 

dr 


M 


(3-2) 


N  -   1/k  '-f-     N 


1/kr  f^f  N 


-1/kr  3-(r|S) 

dr  dr 


-1/r 


,d!g 


(3-3) 


II'-  r  <•  is  used  as  constant  vector  the  z  unit  vector  a  ,   In  general,  how- 
ever, the  divergence  of  unit  vectors  is  not  zero,  and  thus  unit  vectors 

■>"■    not  fixed,  e  g  ,  in  spherical  coordinates  div  a   -  2/r 

— r 
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A  more    usable    form  of    the   last    is 


Nz  =   1/k    [k2g  +  -£]    , 


(3-4) 


The   scalar    solution   is  well   known: 


f   ~   Z  (/k     -  3  r)  [£   cos   (3z   +   C  sin   (3z]  [A  cosmcp   +  B  sin  mcp] 


where  Z      is   a   linear   combination   of  J    and  N   Bessel    functions. 

m 


Z=b 


Z=0 


r? 


Figure  3=   Parallel  Plate  Regi 


on 


The  separation  constants  are  3  and  m.  The  complete  solution  is  a  double 
summation  over  the  possible  values  of  (B  and  m„  The  two  vector  solutions 
give  rise  to  two  types  of  fields;  TE  and  TVL   For  the  TM  modes,  let 

E  -  1/k-N  ,   H  =  -J-  M  • 

"g     -    IK   g 


while  for  the  TE  modes 


E  ■  Mf) 


H  =  -  N  . 

-     n  -f 


The  f  and  g  scalars  differ  only  in  the  constants  which  are  as  yet  unspee 
ified.    The  boundary  conditions  require  that  the  tangential  electric 
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field  be    zero   at    the    two  metal   plates.      So    for    z^0,bE(p^0=E  For 


the  TM  case    this  means   that 


(£  cos   |3z  -  £   sin   (3z)  =  0 

|z   =  0,b 


From  this  it  is  seen  that 


a 

C  =   0,   sin  3b  -  0,   and  3  -  ~,    I   =  0,1,2,3,   (3-5) 

b 

The  other  constants  are  not  determined  yet  so  B,   may  be  put  equal  to  1  with- 
out any  restrictive  effect.   The  TE  case  is  analogous,  giving 

(£  cos  Bz  -  C   sin  3z) |        =  0 

z  -  0,b 

and  £,   =  0,  C  -    1,  and  6  is  the  same  as  above. 

These  values,  together  with  the  equations  (3-3,4)  give  the  field  ex- 
pressions for  the  permissible  modes  in  the  parallel  plate  line. 
In  full  form  these  expressions  are 


TM. 


8    3 
E_   =  — —  —  Z_   sin   3z    [A  cos  m<p   +  B  sin  mcp]  (3-6) 

k2    dr  ^ 

^m   "   ~,9"  2w.   sin   3z    [B  cos  mcp  -  A  sin  mcp]  (3-7) 

v  k    r     '" 


k2  :3: 


Ez — —  Zm  cos   3z    [A  cos  mqp  +  B  sin  mcp]  (3-8) 


H_  -*—  Zm  cos  3z    [B  cos  mcp  -  A  sin  mcp]  (3-9) 

1  T^kr     m 

II  /-      cos    3z    [A  cos  mcp   +   B    sin   mcp]  '    (3-10) 

^  nk     )r 


TE: 

Er   =  jr  Zjj,  sin   (3z    [D  cos  mcp  -  C  sin  mqp] 

E^  =  -  —  Zm  sin   |3z    [C  cos  mqp   +  D  sin  mqp] 

j(3       Q 

H_  =  —  —  Z^  cos    (3z    [C  cos   mcp  +  D  sin  mqp] 
1        nk  or     m 


«, 


jm0 
nkr 


Zj^  cos   (3z    [D  cos   mcp  -  C  sin  mqp] 


•   /  i    2  q  2  \ 

H     =   : Zm  sin   (3z    [C  cos  mqp   +  D  sin  mcp] 

L  T)k      m 
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(3-11) 
(3-12) 

(3-13) 

(3-14) 

(3-15) 


The  TE,  TM  nomenclature  here  refers  to  transverse  fields  with  respect  to 
the  z  axis.   Neither  set  is  transverse  to  the  r  axis. 

3,3  Thin  Radial  Lines 

For  both  TE  and  TM  modes,  the  cutoff  wavelength  is  given  byA*«j),(24) 


K   =  == 


1 

\lV  +  |  m  1 
2b           \2kt 

2 

% 

(3-16) 


After  rearrangement,  the  condition  for  propagation  is 


r 


il2  + 


471 


4/b1: 


J4 


(3-17) 


If  'i  i   0  and  for  m  -  0,  propagation  occurs  for  b/X  >  /C/2;  this  criterion 
is  the  same  as  that  for  a  plane  wave  between  parallel  planes.   If  i>   =  0 
there  is  propagation  for  2n;r/X  >  m,  and  it  may  be  noted  that  this  con- 
dition involves  the  radius;  a  wave  may  be  evanescent  close  to  the  source 
and  propagating  farther  out.   These  modes  vary  with  the  azimuthal  angle. 
For  the  assumed  distribution  of  tangential  E  in  the  gap,  it  can  be 
seen  from  Eqs.  3-6  to  3-15  that  (3  =  0;  the  equations  reduce  to  a  simple 
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triplet.  The  transverse  fields  (transverse  to  the  radius)  can  be  written 
in  terms  of  mode  voltages: 

Ez  =  2  Vm  (kr)  $(<p)  (3  18) 

z   m  "' 

KL  =  2  -f  Vm(kr)  $(<p)  (3-19) 

For   the   even  modes, 

Vm(kr)    =   An  ^    (kr),      $m   =    cos   m<p  ,  (3   20) 

while    for    the  odd  modes,, 

Vm(kr)    =   Bm  Z^   (kr),      $m   ■   sin  mcp  ,  (3-21) 

Next  the  internal  (radial  line)  field  will  be  expressed  in  terms  of 
the  assumed  distribution.  Equation  3-18  is  to  be  used  for  Ez  (or  Eq)  in 
Eq.  2-28  and  the  integration  performed.  It  is  possible  to  integrate  the 
step  function  by  assuming  the  gap  width  is  quite  narrow,  sin  9-1. 
Rather  than  perform  this  evaluation,  a  more  elegant  scheme  involving  the 
Delta  function  will  be  used;  this  allows  the  effect  of  the  gap  spacing  to 
be  observed 

3  4  Evaluation  of  the  Integrals 

/on 
Infeldv^°'  has  shown  that  the  use  of  a  constant  gap  distribution  is 

equivalent  to  using  the  Dirac  Delta  function^"'  and  evaluating  the  field 

at  the  edge  of  the  gap   He  treated  a  spherical  antenna  with  symmetrical 

excitation  and  found  the  input  admittance  at  the  center  of  the  gap  for  a 

t * p  function  voltage   This  was  in  the  form  of  an  infinite  series  which 

then  truncated  for  computation   He  also  solved  for  the  admittance 
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using  a  Delta  function,  evaluating  the  field  at  the  gap  edge.   If  the 
field  were  evaluated  at  the  center  of  the  gap,  the  series  would  diverge 
to  infinity  which  is  what  one  should  expect  from  a  zero  thickness  genera 
tor.   The  Delta  function  voltage  was  integrated  by  considering  a  func- 
tion of  a  dummy  variable  s  which  approached  Delta  as  s  -*  °°.   The  func- 
tion had  an  integral  of  one,  and  was  continuous  and  dif ferentiable . 
This  was  integrated  and  the  limit  taken.   The  interesting  result  was 
that  the  admittance  series  obtained  had  the  same  value  and  the  same 
truncation  error  as  that  obtained  from  the  step  function. 

Put 

Ee  -   Ez  6(^  -  6)  (3-22) 

where  the  line  thickness  b  has  been  made  equal  to  1    Equations  18  and 
22  of  this  chapter  are  inserted  into  Eq.  2-28  giving 


V  =    ema     - J2V J      aE  Vi*i6(|-e)a0  -e  sinGdedq).   (3-23) 
nm   2rcn(n  +  1)  °   „aD  i   1  x  2    -9  -nm 


The  Delta  function  reduces  the  integration  to  the  value  of  the  integrand 
at  9  -  |;  trigonometric  function  orthogonality  condenses  the  sum  on  i  to 


one  term,  with  the  result 


v  (ka)  .  ^o-Liflj-W!^  (3_24) 


nm         n(n  +  1)   *m  (<p) 


where  the  a   •  e   and  the  $m  are  both  even,  or  both  odd, 
-9   -nm         m 

Equation  2-10,  i.e.,  the  9  component,  is  simply  the  Fourier -Legendre 
expansion  of  the  radial  line  field  (strictly  speaking,  a  Fourier-deriva- 
tive of  Legendre  expansion). 
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For  a  given  primary  feed  configuration,  the  gap  voltage  as  a  func- 
tion of  angle  depends  on  the  antenna  impedance  per  unit  width  of  periph 
ery,  which  is  also  a  function  of  angle,  while,  on  the  other  hand,  the 
antenna  impedance  is  a  function  of  the  gap  driving  voltage.   A  better 
understanding  of  why  the  problem  cannot  be  separated  into  two  distinct 
parts:   1)  finding  the  gap  voltage  for  a  specified  radial  line  geometry, 
and  2)  writing  the  antenna  fields  for  the  voltage  found  in  part  one,  may 
be  obtained  by  following  the  simple  example  of  a  two  element  array  with 
each  element  fed  by  a  transmission  line. 

Take  an  array  of  two  monopoles  of  fixed  length  and  spacing,  each 
excited  by  a  transmission  line  from  a  source.   (See  Fig.  4.)   If  the 
base  currents  Ii,  Io  are  known,  the  pattern  of  the  array  can  be  found 
(analysis)  and  if  a  desired  pattern  is  realizable,  the  currents  to  pro- 
duce it  can  be  found  (synthesis).   The  thesis  here  is  that,  disregard- 
ing the  process  of  relating  the  base  currents  to  the  radiated  field 
and  vice  versa,  the  analysis  problem,  including  the  transmission  lines, 
is  considerably  more  difficult  than  the  synthesis   The  latter  simply 
involves  choosing  a  transmission  line  to  produce  Vi  across  the  known 
impedance  Z^.   For  the  analysis  case  given  the  lengths  -ti  and  ^9  anc' 
characteristic  impedances  Zqt>  Zqo  of  the  lines,  the  source  voltage, 
and  the  self  and  mutual  impedances  of  the  antennas  (which  are,  of 
course,  independent  of  the  currents),  find  the  base  currents. 
The  input  impedance  of  each  element  depends  on  the  current  flowing 
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k 


Figure  4       Two  Monopole  Array 


in  the  other  element 


Vi        li   LXx    +   I2   ^12 


V2        li   Z12   +   1 2  Z . 


(3=25) 


The   base  voltages   and   currents   obtained    from   the    transmission   line 
equations    are   given   by 

V  =   Vi   cosh   y^i   ~~   li   Z01  sinh  y^i 

V  =   V2   cosh   y^2   ~~   I2   Z0  2  sinh  y^2,  (3-26) 
Eliminating   voltages   to   get    an  equation   in   the  currents 

V  -    [Ix   Zn   +   I2  Z12]    cosh  y^i   -  li   Z01  sinh  y^i 

V  =    [Ix   Zt,2   +   I2  Z22]    cosh  y^2  ~~   I2  Z02  sinh  y^2 

(3-27) 
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These  equations  are  solved  simultaneously  for  the  base  currents  which 
are  used  to  calculate  the  radiation  pattern.   The  solution  of  Eq.  3-7 
is  straightforward  but  tedious.   In  summary,  the  impedance  of  each 
element  must  be  equated  to  the  impedance  seen  into  the  line,  and  the 
equations  then  solved  simultaneously. 

The  next  chapter  will  be  concerned  with  matching  the  antenna 
impedance  to  the  line  impedance,  or,  what  is  equivalent,  matching  the 
magnetic  fields  at  the  aperture  or  gap. 


BOUNDARY  CONDITIONS 


I.I  The  Transverse  Magnetic  Field 

The  boundary  conditions  at  the  gap  require  that  the  transverse  magnetic 
field  be  continuous  across  the  aperture   With  the  assumed  distribution  of 
tangential  E,  the  radial  line  transverse  magnetic  field  has  only  a  cp  com 
ponent.   Hence,  the  boundary  situation  will  be  approximated  by  matching  this 
component  to  the  cp  component  of  the  antenna  field   From  Eq ,  2-11,  the  cp 
component  of  external  field  is 

J 


H, 


V(ka>  a0  *  £nm(0(P) 


(41) 


Inserting  the  values  for  the  mode  functions  and  voltages  gives 


cd    n 

(UJTz: 


Hcp  =  -J/^  /_  — m  [\a   cos  m(p  "  B-  sin  m^ 


n=l  m=0  n(n  +  X) 


ni 


ka  Ka   PAm(Q)  PA"1^)    rn2  (rhj;  Pnm(0)  Pnm(A) 

ka  h, 


(rhn)a 


A  -  cos  6  at  the  gap  edge. 


na 


(4-2) 


Modes  here  involve  both  n  and  m  but  the  radial  line  field  has  only  m  modes 
Consequently,  it  is  desirable  to  reorganize  the  double  sum  with  the  m  sum 
outermost  so  that  a  term-by-term  comparison  could  be  made  with  the  single 
series  of  the  line.   The  sum  is  detailed  in  Fig   5 
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24 


H, 


nm 


n=l  m= 


0   12   3   4 
Figure  5  Sum  By  Rows 

This  summation  is  equivalent  to  the  following  summation  by  columns.  Note 
that  this  does  not  constitute  a  rearrangement,  as  the  elements  in  any  one 
row  (or  column)  have  not  been  rearranged. 


00      03 


Hn 
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unm   '-'oo 


m-U  n=m 
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Figure  6,  Sum  By  Columns 

To  signify  that  the  n  =  0,  m  =  0  term  is  not  included,  the  notation  used 
is 


00      00 


«• 


nm 


m-0   n=m 


where  n  starts  at  1  if  m  =  0. 
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The   magnetic    field   is    then 


CO 


00 


Hq,  -  -j/n     ZmlAn  cos  m(P  +  Bm  sin  m^ 


(Unm) 


-D 


L   _*  n(n  +  1) 
n  =  m 


ka  hna  P^CQ)  p^(A)  _  m2(rhn);  Pnm(0)  Pnm(A) 

ka  h 


(rhn)a 


na 


(4  3) 


The  inner  series  contains  two  terms;  the  summation  of  the  first  will  be 

designated  S(m,a),  the  second,  T(m,a)0   The  value  of  these  series  depends 

on  a  and  m  alone.   Examination  of  the  convergence  and  values  of  these 

will  occur  later.   Equation  4-3  becomes 

oo 

Hcp   =   ~JM   /       Zm    (S  ~   T)     f\i   cos   m(P   +   Bm   sin   m^     • 
m=D 

Equating   this   to   the   radial    line    field   (3-18)    yields 


oo 


3 


0  =  )   [Zm  (S  -  T)  -1/k  ;£-  ZJ  [A,,,  cos  mcp  +  F^  sin  mcp]    (4-4) 
nFD 

which  must  hold  for  arbitrary  radii  a  and  coefficients  Am,    Bm.   This  is 
possible  only  if  each  term  vanishes.   Thus 


k  Zm  [S(m.a)  -  T(m,a)]  =  §-   Z 


(4-5) 


for  each  m  =  0,  1,  2,  ...   It  will  be  recalled  that  Zm  is  a  linear 
combination  of  Bessel  functions;  since  the  A,  B  coefficients  are  yet 
unassigned,  only  one  constant  need  be  used  in  7^,    the  ratio  is  the 
quantity  of  importance.   Let 

Zm(ka)  -  wm   Jm(ka)  +  Nm(ka).  (4-6) 
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1.2  The  Impedance  Matching  Function 

Equation  4-5  allows  the  TJm  coefficients  to  be  found  for  given  a  and 
m.   It  may  be  recalled  from  Eq.  4-3  that  the  S(m,a),  T(m,a)  are,  in 
general,  complex.   Calling  the  real  and  imaginary  parts  of  (S  -  T)  by 
S  ,  S  and  the  real  and  imaginary  parts  of  TJm  by  HT   ,  w  : 

S1  =  Re  [S(m,a>  -T(m,a)]     S2  =   Im  [S(m,a)  -T(m,a)]      (4-7) 

TS'   =  Re  TJm     ^  -  Im  -0^  .  (4-8) 

The  solutions  for  -®~    are 


r»   _       (SX-  Jm)(SX-Nm)   MS2)2  JmNn 

(S1^- J^)1  +  (s2Jm)2 


(4-9) 


2 


S2 [J:Nm  -  JmN'] 


t^'   =  — -i — "  m  m  a     m     2 2-  ■  (4-10) 

m        (s1^- Ji)"  ^  (S"Jn)a 

It   may   be   easier   to    find  tsl      in    terms   of  THL    . 

'  ,  m  m 

yi.gl\-t  -V  s%  ■  (4-ii) 

The  quantity  of  main  interest  is  Z  which  turns  out  to  be 

Z  (ka)  -  [J"(Sl  "  jS2)  "  Jm^JmNm  ~  JmNm)  {^u) 

(s1^  -  Jm)2  ♦  (s2Jm)2 

After    replacing   the  Wronskian   by   2/rcka,    the    final    form   is   obtained. 

Zm   ■■   -2-       VSl   "J   S2)   -  Jm       .  (4.13) 

nka  (s1^- Jm)a  +  (s2Jm)2 


Primed   Bessel    functions    indicate   derivatives  with    respect   to   argume 


nt.  I 
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The  Z   function  is  related  to  the  relative  proportion  of  radiation  from  one 
m  mode,  assuming  that  all  modes  are  excited  in  equal  amounts.   The  actual 
amount  of  each  mode  excited  is  dependent  on  the  feed  arrangement;  i.e.,  on 
the  A,  B's. 

It  is  next  necessary  to  arrive  at  some  numerical  values  for  the  S(m,a) 


and  T(m,a)  series  which  are  repeated  here, 

oo 


S(m>a)  .  V       (2n  +  l)(n-m)!  ka  hna  Pn"(0)  P^ 
n^,*        2n(n  +  l)(n  +  m) !  (rhn)a 


(A) 


(4-14) 


00 


T(m,a) 


I 


n=m 


(2n  +  l)(n  -  m)!m2  (rhj'  Pnm(0)  Pnm(A) 


n"a   n 


2n(n  +  l)(n  +  m)!  ka  hn; 


*  (4-15) 


The  rapidity  of  convergence  may  be  examined  through  the  use  of  asymptotic 
formulas  for  large  n.   From  Stirling's  formula^'', 

2m 


(2n  +  l)(n-m)! 


for  large  n. 


(4-16) 


2n(n  +  l)(n  +  m) !    n2m  +  1 

From  Morse  and  Feshbach'™) >  the  asymptotic  form  of  the  Hankel  function 
for  large  order  may  be  found. 

I^(2)(x) 


_2_    f2jl]n    (e)-n   +  &c   +  Y2KJ 
>l7in      x 


(4-17) 


and  hn(2)(x)    =   [^Hn?i/2   (x)    so   that 
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h   (2)    (x)   .  1   lTr2n]n    (e)-n   +  fee   +  ft*  (4_lg) 

x  nI  e    [  x  J 


Thus 


ka  h 


na 


h(2) 
na 


<rhn>a 


ka   h 


(2) 

n  -   1,  a 


-  n 


(2) 
4ia 


ka 
n 


(4-19) 
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The  asymptotic  formulas  for  the  associated  Legendre  functions  for  large 
lues  of  index  n  are,  from  Magnus  and  Oberhettinger^  ^' : 


va 


(-l)m+1 


2-  nm+Y2  sin[(n    +  V2)    9   +   S&  -  H] 
S-  Pnm(cos   0) i5 2 4_  (4.20) 


36     n  e^  sin^  9 

(-l)m    f2~nm-l/2  cos[(n   +  %)    9  +  ^  -  ?] 
>J  tc                                              2        4 
Pnm(cos   6) (4-21) 


e 


"^  sin^  9 


These  approximations  are  for  n  —  °°  and  for  9  near  n/2.      Insertion  of  the 
foregoing  asymptotic  forms  into  Eqs .  4-14  and  4-15  for  the  series  shows 
that  the  ntn  term  of  both  varies  as  1/n  plus  a  slow  sign  change;  the 
latter  produces  convergence.   These  series  depend  on  both  a  and  m  and 
therefore  many  computations  of  tedious  nature  would  be  needed.   An 
algorithm  is  needed  to  ease  this  problem.   An  attempt  was  made  to  find  a 
closed  form  for  the  series.   The  general  term,  involving  the  product  of 
two  associated  Legendre  functions  bears  resemblance  to  the  Addition 
Theorem  for  those  functions.   The  latter,  however,  contains  a  summation 
on  m  whereas  the  series  at  hand  involves  n.   Another  possibility  sug- 
gested was  the  method  of  replacing  the  sum  by  a  contour  integral  where 
each  term  is  the  residue  of  the  function  at  a  pole  on  the  real  axis. 
The  trick  is  to  replace  this  contour  by  one  including  the  other  singu- 
larities of  the  function  than  the  poles  on  the  axis.   Often  this  pro- 
duces a  much  more  rapidly  converging  series.   This  technique  is  elabo- 
rated i  ii  Appendix  V  ■ 
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4  3  Increasing  the  Rapidity  of  Convergence 

Another  technique  that  is  often  useful  in  situations  of  this  type  is 
that  of  adding  and  subtracting  a  suitable  convergent  series  to  the  series 
under  study.   The  new  series  is  added  term  by  term,  and  subtracted  as  a 
sum.   This  new  series  is  chosen  to  make  the  resulting  series  converge 
more  quickly.   To  this  end  two  series  are  defined. 


oo 


m 


n=m 


(2n  +  l)(n  -  m)!  Pnm(0)  Pnm(A) 
2n2(n  +  l)(n  +  m)! 


00 


(4-22) 


n=m 


00  00 

y  (2n  +  i)(n  i  m)!  pnm(o)  pnm(A)  .  y 

L^~*  2(n  +  l)(n  +  m)'         *—=- * 

n-m         ^v"  L>\"        »!/•  n=m 


(4-23) 


Subtracting  and  adding  these  to  Eqs.  4-14  and  4-15,  respectively,  produces 


S(m, a)  =  -ka  J  S 


m 


/  nhna   +  , 


(4-24) 


T(m,a)  = 


(rhn>a  ,  - 
nhna 


(4-25) 


From  Eq.  4-19  it  can  be  seen  that  as  n  increases,  the  factor 

nhr 


*na 


(r<Va 


+  1 


decreases  quite  rapidly.   The  original  problem  has,  therefore,  been  split 
into  two  problems:   1)  calculation  of  the  slowly  convergent  Sm  and  Tm 
series  which  are  independent  of  a;  2)  using  the  an,  bn  already  found  in 
the  rapidly  converging  series  of  Eqs.  4-24  and  4-25.   The  rapidity  of 
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convergence  of  the  latter  is  a  function  of  the  argument  ka,  being  more 
rapid  for  small  values,  but  it  is  of  the  order  of  1/n  so  that  only  a  few 
terms  will  be  needed. 

The  calculation  of  the  Sm,  Tm  series  was  made  on  the  ILLIAC,  the 
University  of  Illinois  digital  computer.   This  machine  is  a  high  speed 
single  address  parallel  machine  with  a  fixed  binary  point  and  has  a  high 
speed  electrostatic  memory  storing  1024  forty  bit  words  which  gives  a 
single  precision  of  twelve  decimals.   The  add  time  is  75  microseconds. 
Automatic  conversion  codes  are  available  in  addition  to  a  large  library 
of  commonly  used  subroutines. 

The  calculation  is  described  in  detail  in  Appendix  I;  a  brief  order 
code  for  the  ILLIAC  is  included  in  Appendix  II  which  will  allow  programmers 
to  get  the  gist  of  the  codes.   Results  of  the  calculations  are  given  in 
the  first  appendix.   The  next  step  is  the  computation  of  the  S(m,a), 
T(m, a )  series. 
4  4  Calculation  of  the  Series 

The  coefficients  an,  bn  which  depend  on  n  and  m  have  already  been 

calculated  by  the  computer  codes  of  Appendix  I   The  remainder  of  the 

expression  of  4-24  varies  with  the  antenna  radius  a  and  n.   To  calculate 

the  sum  for  a  given  a  and  m,  the  ILLIAC  was  again  utilized  with  the 

pertinent  details  relegated  to  Appendix  III   Because  no  codes  for  calcu- 

i-i'irit.'  spherical  Bessel  functions  were  available  and  due  to  shortage  of 

tine  for  writing  such  a  code,  values  were  picked  from  tables^    , 

punrhcf]  on  tape,  and  read  into  the  machine.   The  values  of  ka  used  are 
List  ed  i  n  Tab]  e  I 
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ka       a/A  ka       a/A 


0,16 

1/40 

0.31 

1/20 

0.47 

3/40 

0,63 

1/10 

0.78 

1/8 

1.05 

1/6 

1.25 

1/5 

1.57 

1/4 

1.88 

3/10 

2.05 

0.32 

2.20 

7/20 

2.40 

0  38 

2.52 

2/5 

2.83 

9/20 

3.14 

1/2 

3.50 

5/6 

3.76 

3/5 

4.71 

3/4 

5.34 

17/20 

5.80 

0  92 

6.28 

1,00 

TabVe  i.     Values  of  Argument 

The  an,  bn  values  were  also  put  on  tape  to  be  read  in.   The  general 
scheme  of  the  calculation  is  as  follows.   The  master  tape  is  read  in, 
then  the  an  tape.   Then  a  set  of  spherical  Bessel  functions  for  one 
argument  (on  tape)  is  read  in  and  the  sum  found  and  printed.   Addi- 
tional sets  of  SBF  may  then  be  read  in  and  the  sums  for  these  argu- 
ments found.   Available  tables  limit  the  number  of  terms  to  ten  but 
this  gives  ample  accuracy  due  to  the  rapid  convergence  of  the  series. 
Values  of  S(m,a)  and  T(m,a)  were  hand  calculated  (on  a  desk  calcu- 
lator, of  course)  from  these  data  (See  Tables  2  through  5.)   The  real 
and  imaginary  parts  of  the  first  few  have  been  plotted  in  Figs.  7,  8, 
9,  and  10. 

Next  the  S"  and  S  of  Eq.  4-7  were  hand  calculated,  paving  the 
way  for  the  calculation  of  Z^,  Eq.  4-13,  the  goal  of  the  operation. 
This  last  calculation  was  programmed  for  the  ILLIAC,  the  floating 
decimal  subroutine  (Cf.  AppendixIH)  is  ideally  suited  for  such  cases 
where  the  numbers  have  a  wide  range  of  values.   Details  are  found  in 
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Appendix  IV  the  results  are  plotted  in  Figs.  11  and  12   For  smaller  ka, 
the  magnitude  of  Z^  is  largely  controlled  by  T(m,a)  so,  to  preserve  the 
corresponding  phase,  only  T(m,a)  was  used  for  ka  £  2.83   Above  this 
value,  both  T(m,a)  and  S(m,a)  were  used.   This  method  was  found  to  give 
good  results. 

ka 

m  16  .31      .47       63      78     1.05 

0  -  .2636  -  5273  -  8403  -1.179   -1.462   -1.655 

1  -  .1692  -  .3292  -  5032 

2  -  .1319  -  .2563  -  3896 

3  -  .1099  -  2131  -  .3253 

4  -  0943  -  1829  -  .2775 

5  -  0824  -  1597  -  2423 

6  -  0728  -  .1411  -  2139 

7  -  0647  -  1255  -  1902 

8  -  .0579  -  1121  -  1700 

9  -  .0519  -  .1005  -  1524 
10  -  0466  -  .0902  -  1368 


6827  - 

8591 

-1.206 

5242  - 

6520 

-  8884 

4344  - 

5391 

-  7296 

3720  - 

4616 

-  .6234 

3239  - 

4027 

-  .5432 

2869  - 

3554 

-  .4792 

2551  - 

3160 

-  4257 

2279  - 

2823 

-  .3803 

2043  - 

2529 

-  .3407 

1834  - 

2271 

-  .3058 

ka 


m 

1.25 

1.57 

1.88 

2.05 

2.20 

2.40 

0 

-1.642 

-1.737 

-1.993 

-2.182 

-2  366 

-2.616 

1 

-1.492 

-1.965 

-2.254 

-2,304 

-2.324 

-2.372 

2 

-1.071 

-1.384 

-1.723 

-1.927 

-2.114 

-2.361 

3 

-  8732 

-1.  1090 

-1.348 

-1.486 

-1.612 

-1.789 

4 

-  7444 

-  .9407 

-1.135 

-1  244 

-1.342 

-1.476 

5 

-  6479 

-  8169 

-  9829 

-1.075 

-1.157 

-1.268 

6 

-  5713 

-  7194 

8644 

-  9444 

-1.016 

-1.111 

7 

-  .5073 

6384 

7661 

8366 

-  .8989 

9826 

8 

-  4531 

-  5699 

-  6836 

7462 

.8016 

-  .8758 

9 

-  4057 

5101 

-  6114 

-  6671 

-  7164 

-  .7822 

10 

-  .3642 

-  4578 

5486 

-  .5985 

6426 

-  .7015 

Table  2   Real  Part  of  S(m,a)  Series 
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ka 


m 

2.52 

2.83 

3.14 

3.50 

3.76 

4.71 

0 

-2.752 

-2.977 

-3.030 

-3.114 

-3.265 

-3.984 

1 

-2.424 

-2.643 

-2.952 

-3.310 

-3.465 

-3.728 

2 

-2.495 

-2.730 

-2.810 

-2.945 

-3.111 

-3.842 

3 

-1.901 

-2.212 

-2.549 

-2.912 

-3.084 

-3.468 

4 

-1.558 

-1.779 

-2.016 

-2.321 

-2.563 

-3.370 

5 

-1,335 

-1.523 

-1.696 

-1.921 

-2.097 

-2.834 

6 

-1.169 

-1.321 

-1.475 

-1.660 

-1.798 

-2.350 

7 

-1.033 

-1.164 

-1.298 

-1.455 

-1.571 

-2.016 

8 

-  .9206 

-1.037 

-1  154 

-1.292 

-1.393 

-1.774 

9 

-  .8218 

-  .9245 

-1.028 

-1.149 

-1.237 

-1.565 

10 

-  .7371 

-  .8289 

-  .9210 

-1.029 

-1.107 

-1.398 

ka 


m 

5.34 

5.80 

6  28 

0 

-4.075 

-4.264 

-4.596 

1 

-4.190 

-4.372 

-4.405 

2 

-3.960 

-4.161 

-4.495 

3 

-3.946 

-4.144 

-4.208 

4 

-3.586 

-3.837 

-4.176 

5 

-3.376 

-3.623 

-3.765 

6 

-2.791 

-3.160 

-3.546 

7 

-2.343 

-2,609 

-2.918 

8 

-2.043 

-2.252 

-2  495 

9 

-1.790 

-1.960 

-2.141 

10 

-1.595 

-1.743 

-1.907 

Table  2   Real  Part  of  S(mva)  Series  (Cont  ) 


ka 


ka 


ka 


2.52 


2.83 


3.14 


3  50 


3  76 


ka 
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m 

.16 

.31 

.47 

.63 

.78 

1.05 

0 

0005 

,0076 

.0442 

1553 

3644 

.8192 

1 

- 

- 

.  0002 

.0011 

.0040 

.0248 

2 

- 

- 

- 

- 

- 

0004 

3 

- 

- 

- 

- 

- 

- 

m 

1.25 

1.57 

1.88 

2.05 

2.20 

2.40 

0 

.9758 

.9987 

.9854 

1.005 

1.051 

1.169 

1 

.0733 

.2858 

.6476 

.8133 

.9021 

.9517 

2 

0015 

.0089 

0374 

.0741 

.1285 

.2451 

3 

- 

0002 

.0012 

.0028 

.0055 

.0127 

4 

- 

- 

- 

.0001 

0002 

.0005 

5 

- 

- 

- 

- 

- 

- 

4,71 


0  1.274   1.603   1.841   1.928   1.952   2.513 


1 

9606 

.9774 

1.057 

1.305 

1.553 

1.936 

2 

3414 

.6350 

.8485 

.9359 

.9684 

1.525 

3 

.0203 

0611 

.1597 

.3887 

.  6022 

.9629 

4 

0009 

0033 

.0104 

0348 

.0755 

.5863 

5 

- 

0001 

0005 

.0021 

0054 

.0897 

6 

- 

- 

- 

.  0001 

0003 

.0080 

7 

- 

- 

- 

- 

- 

.0006 

8 

- 

- 

- 

- 

- 

- 

m   5  34 


5.80 


6  28 


0 

2  853 

2 

971 

3  183 

1 

2  216 

2 

554 

2  818 

2 

1  858 

1 

997 

2  194 

3 

1  249 

1 

576 

1  827 

4 

8828 

9958 

1.233 

5 

3582 

6576 

8613 

6 

0459 

1380 

3520 

7 

0044 

0155 

0513 

8 

0003 

0015 

0056 

9 

- 

0001 

0005 

L0 

- 

- 

- 

Table  3. 

Imaginary  Par 
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ka 

m 

16 

31 

47 

63 

78 

1.05 

1 

-  7 . 040 

-  3  544 

2  266 

-  1  624 

-  1  263 

-   8781 

2 

-  24.94 

-  12.81 

-  8  159 

-  6.180 

-  4  916 

3,522 

3 

-  50.75 

-  26.14 

-  17.01 

-12.77 

-10,25 

-  7  507 

4 

-  82.80 

-  42.69 

-  27  95 

-20.92 

-16  84 

-12,42 

5 

-119  9 

-  61.86 

-  40  61 

-30.35 

-24 , 47 

-18.09 

6 

-161  1 

-  83.12 

-  54  65 

-40.83 

-32  93 

-24.39 

7 

-205.7 

-106.2 

69  88 

-52  16 

-42  13 

-31  19 

8 

-253,0 

-130.6 

85.98 

-64  17 

-51  80 

-38.41 

9 

-302.4 

-156=0 

-102.7 

-76.71 

-61  89 

-45.93 

10 

-353.3 

-182.3 

-120  1 

-89  66 

-72  38 

-53.72 

m 

1.25 

1  57 

ka 
1,88 

2,05 

2.20 

2.40 

1 

7055 

.5268 

4155 

-   3698 

-  .3358 

-  .2978 

2 

-  2.861 

-  2.149 

1  695 

-  1.508 

-  1  369 

-  1.214 

3 

-  6  216 

-  4  801 

-  3  859 

-  3  456 

-  3.149 

-  2.799 

4 

-  10.36 

-  8.126 

-  6  659 

-  6.029 

-  5.547 

-  4.988 

5 

-  15.13 

-  11.94 

-  9  862 

-  8.982 

-  8.304 

-  7.536 

6 

-  20.42 

-  16.16 

-  13  40 

-12  23 

-11  34 

-10  70 

7 

-  26  14 

-  20.72 

-  17  21 

-15  73 

-14  61 

-13.33 

8 

-  32.21 

-  25  56 

-  21.53 

-19  45 

-18  08 

-16  52 

9 

-  38  54 

-  30  60 

-  25  47 

-23.32 

-21,69 

-19.82 

10 

-  45.08 

-  35.82 

-  29  85 

-27.34 

-25.43 

-23.27 

m 

2.52 

2.83 

ka 
3  14 

3  50 

3  76 

4  71 

1 

.2782 

.2365 

2034 

-  .1736 

-   1561 

-  .1108 

2 

-  1.134 

9613 

8277 

7053 

-   6337 

4491 

3 

-  2.616 

-  2.219 

-  1  908 

1  624 

1  457 

-  1.027 

4 

-  4.690 

-  4.027 

3  485 

-  2  974 

2  670 

-  1  879 

5 

-  7  123 

-  6.209 

-  5  453 

4  719 

■  4  265 

-  3  022 

6 

-  9.799 

-  8  610 

-  7.659 

-  6  713 

-  6  131 

-  4.477 

7 

-  12.65 

-  11.17 

-  9.959 

-  8  812 

-  8  103 

-  6  100 

8 

-  15  70 

-  13.89 

-  12  43 

-11.04 

-10.19 

-  7  387 

9 

-  18.85 

-  16  70 

-  14.97 

-13.32 

-12.33 

9  574 

10 

-  22  13 

-  19.64 

-  17  63 

-15  73 

-14  58 

-11  42 

Table  4.   Real  Part  of  T(m,a)  Series 
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ka 

m 

5.34 

5.80 

6.28 

1 

-  .0910 

-  .0792 

-  .0700 

2 

-  .3689 

-  .3244 

-  .2848 

3 

-  .8403 

-  .7257 

-  .6431 

4 

-1.535 

-1.348 

-1.177 

5 

-2.463 

-2.145 

-1.869 

6 

-3.675 

-3.205 

-2.805 

7 

-5.095 

-4.475 

-3.910 

8 

-6.676 

-5.946 

-5.288 

9 

-8.239 

-7.421 

-6.767 

10 

-9.907 

-8.989 

-8.173 

Table  1.     Real   Part  of  T(m,a)  Series   (conto ) 
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ka 


m 

,16 

.31 

.47 

.63 

78 

1  05 

1 

-  ,0094 

-  .0329 

-  .0678 

-  .1059 

-  .1418 

-  .1968 

2 

- 

-  .0006 

-  .0031 

-  0095 

-  .0206 

-  .0561 

3 

- 

- 

- 

-  ,0002 

-  0007 

-  .0039 

4 
5 

- 

- 

ka 

- 

- 

-  .0001 

m 

1  25 

1  57 

1  88 

2.05 

2.20 

2.40 

1 

-  .2292 

-  .2688 

-  2970 

-  .3096 

-  .3196 

-  .3315 

2 

-  .0946 

-  .1692 

-  .2442 

-  .2829 

-  .3148 

-  .3540 

3 

-  ,0100 

-  .0317 

-  .0719 

-  ,1025 

-  .1336 

-  .1798 

4 

-  .0004 

-  .0023 

-  .0081 

-  0145 

-  .0229 

-  .0392 

5 

- 

-  .0001 

-  ,0005 

-  .0010 

-  0019 

-  .0041 

6 

7 

- 

- 

ka 

- 

-  .0001 

-  .0003 

m 

2.52 

2.83 

3  14 

3  50 

3.76 

4,71 

1 

-  .3378 

-  .3526 

-  .3647 

-  ,3765 

-  .3839 

-  ,4046 

2 

-  .3756 

-  .4256 

-  .4698 

-  .5137 

-  5415 

-  .6211 

3 

-  .2091 

-  .2866 

-  ,3620 

-  .4429 

-  .4964 

-  .6594 

4 

-  .0522 

-  .0975 

-  .1590 

-  .2443 

.3106 

-  .5440 

5 

-  0061 

-  .0158 

-  .0347 

-  0737 

-  .1138 

-  .3318 

6 

-  .0004 

-  .0014 

-  .0041 

-  .0113 

-  .0225 

-  .1295 

7 

- 

-  .0001 

-  .0003 

-  .0011 

-  .0027 

-  .0297 

8 

- 

- 

- 

-  .0001 

-  .0002 

-  .0043 

9 

- 

- 

- 

- 

- 

-  .0004 

10 

ka 

m 

5.34 

5.80 

6.28 

1 

-  .4148 

-  .4195 

-  .4262 

2 

-  .6606 

-  .6881 

-  .7061 

3 

-  .7434 

-  .7764 

-  .8425 

4 

-  .6779 

-  ,7885 

-  .8438 

5 

-  .4964 

-  .6114 

-  .7233 

6 

-  .2661 

-  .3850 

-  .5134 

7 

-  .0912 

-  .1688 

-  .2759 

8 

-  .0193 

-  .0468 

-  .0996 

9 

-  .0027 

-  .0085 

-  .0233 

10 

-  .0003 

-  .0011 

-  .0038 

Table  5.      Imaginary   Part  of  T(m,a)   Series 
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5.       OFFSET   FEEDS 

5. I     Wire  Feeds 

It  is  the  purpose  here  to  examine  a  simple  wire  feed  for  the  radial 
transmission  line  and  to  determine  the  effects  of  locating  the  feed  eccen- 
trically.  Let  the  feed  be  a  small  wire  of  radius  c  carrying  a  current,  I. 
Since  both  the  wire  diameter  and  length  are  small  compared  with  wavelength; 
one  may  assume  that  the  current  is  uniform  along  the  length  and  around  the 
periphery  of  the  wire.   Furthur,  assume  perfect  conductivity  of  the  wire 
so  that  the  internal  field  is  zero.   Then^^' 


I  =  2nc  K 


(5-1) 


where  the  wire  is  centrally  positioned.   The  general  radial  line  field  was 
given  by  Eq.  3-18. 


oo 


H^  =  -j/n  2_  Zm[Am  cos  mcp  +  Bm  sin  m^  '  (3"20) 

m=0 

^  is  now  a  known  function;  the  A,  B  coefficients  are  determined  by  com- 
paring these  two  equations  for  H™.   It  is  evident  that  for  this  trivial 
case  all  EL  are  zero  and  also  all  Am  except  AQ.   The  latter  becomes 


jr\l  ,     „,    d  Z0 

—  where  Z00 


°   2nc  Z^o         uu   d  ki 


(5-2) 


Again  the  prime  on  the  Z  indicates  differentiation  with  respect  to  argu- 
ment.  From  the  above,  the  radiated  field  may  be  readily  found. 
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5.2  An  Offset  Wire 

Moving  the  wire  of  Section  5.1  off  center  appears  to  be  a  promising 
type  of  feed.   Let  the  wire  be  situated  a  distance,  d,  from  the  center. 
(See  Fig.  13.) 


0  d 

Figure  13,   Addition  Theorem  Coordinates 
The  coordinates  of  the  wire  are  now  (r^qPi).   The  Bessel  function  Addition 
Theorem^  '    is  useful  for  expanding  the  field  in  the  coordinates  (r,qp). 
The  field  of  interest  is  that  at  r  -  a  because  the  antenna  impedance 
boundary  condition  must  be  satisfied  at  the  gap.   In  this  region,  the 
condition  for  convergence  of  the  first  form  of  the  Addition  Theorem, 
namely,  r  >  d  cos  cp  is  always  satisfied  and  the  series  converges.   This 
form  gives  the  field  as  a  sum  of  modes  or  waves  originating  from  the 
center.   Note  that  the  origin  is,  and  must  be,  excluded  from  the  region 
of  analyticity  of  the  function  as  a  spherical  wave  may  come  from  the 
direction  of  the  center  but  cannot  start  at  the  center,  i.e.,  the  Hankel 
functions  have  a  singularity  at  r  -  0.   The  first  form  of  the  Addition 

Theorem  is 

oo 

/.n(krj)  cos  n<pt  =  )    J^(kd)  Z^  Jr  n(kr)  cos  (i  +  n)cp  ,.    (5-3) 
1.-oo 
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Zn  here  means  a  general  Bessel  function,,  either  Jn  or  Nn  or  a  Hankel,  and 
is  not  to  be  confused  with  the  Zm  impedance  function.   Inserting  Eq.  5-2 
into  the  general  expression  for  the  radial  line  field  (3  18)  gives  the 
primary  electric  field  in  the  (r^cPi)  coordinates. 


jnl  H0(2)  (kr)  dH0(2)  (kr) 
with  Hoo  =  " 


2nc  Hoo 


d  ki 


(5-4) 


Expanding  this  field  in  the  (r,cp)  coordinates  yields 


jnl 


oo 


2nc  Hoo 


EiJ^kd)  cos  icp  Hi(2)  (kr) 


i-0 


for  the  primary  field.  The  secondary  field  can  contain  only  J  functions 
because  the  N  functions  have  a  singularity  at  the  origin.  The  secondary 
field  is  then 


oo 


Ez2  =  -JT]  /        eiai  J^(kr)  cos  it 


(5-5) 


i=0 


the  a^  are  as  yet  unknown  constants.   The  total  field,  which  is  the 
superposition  of:   1)  the  primary  field,  i.e.,  the  field  in  an  infinitely 
long  line,  and  2)  the  secondary  or  scattered  field,  becomes: 

Ez  =  jn  Y    ei  cos  1(P<  V1  !>d)  Hi(2)(k^)  -  ai  Ji<kr)l  .   (5-6) 
4—  znc   Hoo 


The  general  radial  line  field  was  (3-18) 

CO 

Ez   =    )       Zm[Am  cos  mcp   +  ^  sin  mcp] 


(3-18) 


m=0 
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and  it  is  desired  to  find  the  values  of  A  and  B  that  will  make  this  iden- 
tical with  Eq.  5-6.  Due  to  the  choice  of  axes,  even  symmetry  is  enjoyed, 
all  Bjp  are  zero.   Equating  like  terms  in  cos  mcp, 


jner 


I  Jm  »«i)   (2)  _ 
2nc  Hoo  ^ 


am^m 


^m  \i   • 


(5-7) 


This  must  be  valid  for  all  r  >  d  so  the  J  and  N  portions  must  independ 
ently  be  zero,  giving 

IemJm(kd) 


\,  =  ^ 


'rrr-'m 


2^c  Hoo 


(5-8) 


x'*  Jm  (kd> 

am  =  (1  +  jtsl)    25 

m       J  m   2nc  Hoo 


(5-9) 


This  result  behaves  about  as  one  might  expect:   as  d  -*  0  the  amounts  of 
higher  modes  excited  diminishes  giving,  in  the  limit  d  =  0  only  the 
dominant  mode  because  all  Jm(0)  are  zero  except  J0. 

The  constant  Hoo  is  generally  complex  but  may  be  simplified  for  the 
small  argument,  c  <<  A. 


-^-Ho<2>(kr) 
d  kr 


■  T-f"  [Jo(kr)  -  j  No(kr)] 
c   d  kr 


Ji(kc)  +  j  IMkc). 

(5-10) 


For  small  arguments, 


Jr   ~  /4  kc,   Ni  ~ 


_2_ 

nkc 


which  gives 


Hoo  "  -%  kc  -  j  -2- 

nkc 


(5-11) 
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Since  kc  <<  1, 

Hoo  1  -   j  -2-  .  (5-12) 

nkc 

Equation  5-8  is  then 

Am  =  J^k  *  em  I  Vkd)  ■  (5-13) 

It  is  seen  that  the  A' s  are  all  imaginary,  a  fact  which  will  simplify  the 
calculations.   For  purposes  of  calculating  patterns,  all  of  the  constants 
can  now  be  lumped  into  a  single  constant  and  equated  to  one.   The  desired 
pattern  is  then  found  by  inserting  Eq.  5-13  into  Eqs.  3-20,  3-21,  and  - 
3-24,  removing  superfluous  constants,  and  then  using  these  values  in  the 
field  equations  (2  10).   Typical  patterns  would  be  for  a  fixed  value  of 
9,  giving  the  electric  field  amplitude  as  a  function  of  <p. 
5.3  Calculation  of  Patterns 

Scrutiny  of  the  field  equations  reveals  that  it  is  necessary  to 
calculate  the  terms,  which  are  complex,  of  a  double  series  and  sum  these 
terms  for  each  value  of  cp  (or  9)  needed,   This  involves  considerable 
tedium  and  is  a  fertile  field  for  human  errors,  so  the  calculation  of 
patterns  was  programmed  for  the  ILLIAC. 

The  functions  needed  are:   cylindrical  Bessel  functions  of  the  first 
kind,  spherical  Hankel  functions,  the  Z^   impedance  functions,  and  the 
associated  Legendre  functions  and  their  derivatives.   The  Bessel  functions 
needed  were  relatively  few  in  number,  hence,  were  punched  on  tape  from 
standard  tables  as  described  in  Appendix  M   A  code  was  written  by  the 
author  for  associated  Legendres  and  derivatives  since  a  large  number  are 
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needed  for  each  pattern.   This  code,  discussed  in  AppendixVI,  calculates 
Pnm  (cos  9)  for  a  given  0  and  n,  storing  values  for  each  m  from  0  to  n  in 
a  preset  set  of  locations  in  the  memory.   A  recursive  technique  is  used; 
the  code  will  also  calculate  the  derivative. 

Two  codes  were  written  for  the  calculation  of  offset  patterns.   One, 
written  by  Mr.  Earl  Schweppe  of  the  Antenna  Laboratory,  is  an  elaborate 
production  which  sums  the  double  series,  and  allows  the  operator  to 
select  and  calculate  any  0  or  cp  pattern  with  any  combination  of  parameters. 
Low  speed  is  the  limitation  of  this  code;  the  slowness  is  due  to  the 
versatility  and  to  summing  the  double  series  for  each  pattern.   The 
second  code,  by  the  author,  breaks  the  calculation  into  two  parts.   The  n 
summation,  which  depends  on  ka  and  m  alone,  is  summed  once  and  for  all  by 
the  synthesis  code  of  the  next  chapter.   These  results  are  read  in  for  a 
desired  ka  and  the  second  code  quickly  calculates  the  pattern.   The  code, 
however,  is  limited  to  qp  patterns  for  which  0  =  90°. 

Some  results  are  presented  for  two  values  of  offset  in  the  next 
pages.   Figures  14,  15,  16  depict  azimuthal  patterns  for  0  =  90°  with  the 
feed  a  distance  from  the  center  of  0.4  times  radius.   Figures  17  through 
20  are  for  an  offset  of  1.0,  i.e.,  the  feed  is  located  at  the  edge.   In 
each  case  the  feed  is  on  the  line  qp  -  0.   These  data  indicate  the  wide 
variety  of  patterns  that  can  be  obtained  from  these  two  variables,  radius 
and  offset.   It  would  take  a  large  number  of  patterns  to  indicate  all  the 
possibilities  even  for  the  range  ka  £  6.28  (diameter  <  2X). 
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ko  =  1.57 


ko  --  1.88 


ko  =  2.20 


ko  =  2.52 


cp  patterns,      9  polarization,      0  =  90° 
Figure   14       Theoretical    Patterns   for  0.4  Offset 
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ko=  2.83 


ko=  3.14 


ka=  3.50 


ko=  3.76 


cp  patterns       9  polarization       9  =    90° 
Figure   15       Theoretical    Patterns  for  0.  i+  Offset 
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ka  =  4.7l 


ka=5.34 


ka=5.80 


ka=6.28 


qp  patterns^      9  polarization,      6  =   90° 
Figure   16       Theoretical    Patterns  for  0  4  Offset 
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kci  =  0.3l 


ka=0.63 


ka=l.05 


ka  =  l.57 


cp  patterns,      G  polarization.    9       90° 
Figure    17       Patterns   for    I    0  Offset 
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ko  =  q  7i 


feM5.34 


ka=  5  80 


ka  =  6.28 


cp  patterns,      0  polarization 

Figure  20   Patterns  for  1.0  Offset 
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Experimental  measurements  were  conducted  using  the  equipment  and  setup 
described  in  Chapter  8.   The  antenna  was  a  hemisphere  of  spun  copper  and 
six  inch  diameter  with  a  brass  plate  soldered  across  the  base,   The  feed 
was  a  standard  BNC  fitting,  with  a  1/16  inch  gap.   Measurements  were  taken 
from  1  to  4  kmc  and  are  the  dashed  curves  of  Figs.  18  and  19,   Several 
factors  limit  the  correspondence  of  these  with  the  theoretical  patterns. 
Several  gap  thicknesses  were  tried  and  it  was  found  that  at  the  higher 
frequencies  the  patterns  varied  appreciably  with  the  gap  dimension. 
Ideal ly;  this  dimension  would  be  as  near  zero  as  possible  but  the  smaller 
the  gap  the  more  critical  the  antenna  is  to  errors  in  spacing  around  the 
antenna  base.   Thus  the  aperture  plate  and  antenna  base  must  be  perfectly 
flat  and  parallel,  a  difficult  condition  to  produce.   Due  to  the  finite 
size  of  the  receiving  horns,  the  boom  angle  0  could  not  be  90°,  the 
closest  value  obtained  was  83°.   This  also  yields  a  slight  divergence. 

From  the  patterns  taken,  it  appears  that  the  gap  introduces  a  fre- 
quency shift  of  about  five  percent.   That  is,  the  theoretical  patterns 
corresponded  most  closely  to  those  measured  at  a  lower  frequency. 

Variation  of  a  typical  pattern  with  the  offset  is  shown  in  Figs.  21 
and  22  where  patterns  are  given  for  offsets  of  0.0,  0,2,  0.4,  0.6,  0.8, 
and  1  0  times  the  radius.   A  smooth  transition  from  the  circular  pattern 
at  the  center  to  the  multi-lobed  pattern  at  the  edge  may  be  seen- 
5.4  Radiation  Resistance 

The  radiation  resistance  for  the  antenna  configuration  of  the  last 
sections  is  readily  found  by  inserting  the  far  field  expressions  into  the 
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kd  =  0 


kd=0.2  ka 


kd=0.4  ka 


kd  =  0.6ka 


cp  patterns,      9  polarization,      9  =   90°,      ka  =  2.83 
Figure  21        Pattern  as  a  Functicr.  of  Offset 
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kd-0.8  ko 


l<d=  ka 


qp  patterns,      6  polarization,      8  =  90°,      ka   *  2.83 
Figure  22       Pattern  as  a  Function  of  Offset 
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integral  of  Poynting  vector,  and  then  evaluating  the  result.   This  gives 

the  radiation  resistance  at  the  wire  feed  point.   In  Chapter  2  the  modes 

of  the  antenna  were  shown  to  be  completely  orthogonal  so  now  each  TE  and 

TM  mode  may  be  taken  separately.   Let  S  be  a  sphere  of  large  radius  r. 

The  radiated  power  for  one  mode  is  given  by 

phm  '  X  /.  E  x  H  dA  =  %  Is   [EqH^  -  E^lR2  sin  0  dBdcp  .    (5-14) 

Equation  5-13  indicates  that  only  \nm   (even)  coefficients  exist;  the 
terms  of  the  integrand  are  then,  for  TM  waves 

EqH,,   -  -^L  Anm  hn    (rhn')(Unm)2[Pnm']2   cos2  mq> 

(5-15) 

E*H9   =    TikimL88  A^   H"    (rhn)(Unm)2[Pn^2    ««'   »*    ■ 
The   power    flow   integral    is    then 

p       .   r[hn    (rhn)   Anm|    ;2^(Unm)2    {[pnm,]2   cos*   m(p  +   _dC [Pnm]  2sin2   m<p} 
ii\i  2r\k  sin   y 

x   sin   0  d0dqp.       (5-16) 
The   integral    is    (2-23);    the   result    from  Chapter  2  may  be  used   to   give 

P 


i  2  I 

nn(n   +   1)    r[hn    (rhn)'    Anml  (5-17) 

TM  T~^ 


m 


for   the   power.      For   the   spherical   Bessel    function    the    large   argument 

forms   are 

h    (2)(kr)   -  1/kr   eJkr   (-j)n+1    , 


(5-18) 


(rhn)'    -  eJkr   (~j)n 
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-hich   allow   (5-17)    to   be    reduced   to 


nn(n   +   1)   A^j 
*B1       ~"  -TT- 


em  ^k' 


(5-19) 


From   (3-24)  and   (5-13), 

=  jem  r,k2a  I  JJkd)  2^  Unm  Pnm(0) 
^m  4n(n   +    1)    (rhj' 


(5-20) 


n'  a 


Th 


e   radiated  power   is    finally    found   to   be   from   the   last    two   equations 


2      2      _2 


PB!     =     ^ 


em7ik a     I'-Ukd)'    [Unm  P:m(0)] 


ii *  n 


16  n(n   +   1) 


ZJka) 


(*KK 


(5-21) 


If   the    radiation   resistance   is   defined   as   R  =   P/I      then 


_        e^Eka  Jm(kd)  UnmPnm(0)] 


16  n(n  +  1) 


<rfcn)a 


(5-22) 


Following  an  exactly  analagous  scheme  for  the  TE  modes  leads  to  the 
radiation  resistance  given  below. 


RTE 


Tt[m  Jm(kd)  Unm  Pnm(0)] 
8  n(n  +  1) 


*na 


(5-23) 


Again  the  calculations  are  tedious  and  lengthy;  again  a  program  was 
written  for  machine  calculation  of  these  values.   Spherical  Bessel 
functions  and  Z   functions  for  a  given  ka  are  read  in,  then  the  Jm  for 
the  kd  are  read  in;  the  operator  has  the  choice  of  printing  each  term  or 
just  the  sum   Additional  ka  and  kd  parameters  are  readily  handled.   Re- 
sults are  given  in  Tables  6  and  7  where  it  may  be  noted  that  the  radiation 

i  itance  goes  through  sharp  maxima  and  that  the  offsetting  of  the  feed 
introduces  more  oscillations  in  the  radiation  resistance  versus  frequency. 
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ka ^        .31      .63     1.25     1.57     1.88 


kd  = 

0  ka 

.0046 

.0204 

.1110 

.2770 

.9383 

.6ka 

.0817 

.1768 

.5266 

l.Oka 

.0044 

.0167 

.0480 

.0682 

.0957 

ka— 

-+ 

2.05 

2.20 

2.40 

2.52 

2.83 

kd  = 

0  ka 

1.331 

.8335 

.3529 

.2353 

.1190 

.6ka 

.5708 

.3678 

.1421 

.0929 

.0977 

l.Oka 

.0709 

.0350 

.0302 

.0341 

.0464 

ka 

— 

3.14 

3.50 

3.76 

4.71 

kd  ■ 

0  ka 

.0831 

.0673 

.0645 

.1895 

.2ka 

.1192 

.9997 

.3925 

.1737 

.4ka 

.1851 

3.135 

.8698 

.  3299 

.6ka 

.2053 

3.409 

.8323 

.5956 

.8ka 

.1481 

1.699 

.2433 

.4627 

l.Oka 

.0657 

.2240 

.0314 

.0912 

ka 

-*- 

5.34 

5.80 

6.28 

kd  = 

0  ka 

6.571 

.4291 

.1282 

.2ka 

3.666 

.2404 

.2299 

.  4ka 

.3370 

.1345 

.7386 

.6ka 

.8090 

.2613 

2.047 

.8ka 

.9931 

.2183 

1.599 

l.Oka 

.0586 

.0717 

.0701 

Table  6.  Radiation  Resistance  Rrad/^  for  TE  Modes 
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ka — 

kd  =   .6ka 
l.Oka 


31 


OOOl 
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0012 


1.25 


0027 
0057 


1.57 


0035 
0073 


1.88 


0044 
0092 


kd  =   .6ka 
l.Oka 


2.05 


0065 
0107 


2.20 


,0076 
,0121 


2.40 

.0114 
.0143 


2.52 


0151 
0156 


2.83 


0263 
0192 


ka- 
kd 


.2ka 
.4ka 
,6ka 
.8ka 
l.Oka 


3.14 

.0178 
.0522 
.0695 
.0538 
.0248 


3.50 


2722 
8977 
9841 
4938 
0685 


3.76 

.0854 
.2127 
.2096 
.0667 
.0138 


4.71 


0165 
1164 
2350 
1791 
0345 


ka- 
kd 


2ka 
.4ka 

6ka 
.8ka 

1  Oka 


5.34 

.0097 
.0478 
,0723 
.0366 
0154 


5.80 


0118 
0416 
0911 
0899 
0305 


6.28 

.0322 

.2907 

.9298 

7180 

0250 


for  kd  =  0,  only  TE  modes  are  excited, 


Table  7   Radiation  Resistance  ^ra(\/r\   for  TM  Modes 
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This  is  quite  different  from  the  resistance  or  conductance  at  the  gap, 
which  exhibits  a  stair-step  hehavior.   The  radial  line  impedance  trans- 
former destroys  this  characteristic  and  substitutes  its  own  sharp 
resonances . 


6   PATTERN  SYNTHESIS  USING  DIELECTRIC  LOADING 

6  I  Continuous  Non-Homogeneity, 

Offset  feeds,  treated  in  the  previous  chapter  can  offer  only  a  limited 
panoply  of  patterns   A  method  of  wider  applicability  involves  loading  the 
radial  line  with  a  material  whose  dielectric  constant  and/or  permeability 
varies  with  the  azimuthal  angle,  <p    By  this  means,  a  variety  of  patterns 
might  be  produced   A  variation  on  this  theme  would  use  a  ferrite  slab 
filling  the  radial  line,  and  would  utilize  a  magnetic  field  produced  by  a 
ring  of  electromagnets  around  the  radial  line  to  change  the  parameters 
of  the  ferrite  material   Then,  by  electrically  rotating  or  sweeping  the 
applied  magnetic  field,  the  antenna  pattern  would  be  correspondingly 
rotated  in  azimuth   This  may  have  military  applications  in  DF  or  search 
where  electrical  scanning  is  desirable 

Since  the  ferrite  loading  is  similar,  analytically,  to  dielectric 
loading  only  the  latter  will  be  examined;  the  results  can  be  immediately 
extended  to  other  cases    It  is  meant  here  that  the  ferrite  materials 
would  be  used  in  that  part  of  the  frequency  range  where  the  gyrator 
effect  is  small   the  ferrite  would  offer  an  enhancement  of  field  by 
U,  with  u.  being  a  function  of  the  applied  static  magnetic  field 

Consider  a  radial  line  filled  with  lossless  dielectric  with 
dielectric  constant  e  -  e  (qp)    In  attempting  an  exact  solution, 
I A"  basic  difficulties  are  met  with  and  these  will  each  be  scrutinized 
I ir  t   when  Maxwell's  equations  are  solved  to  eliminate  U,  there  results 
the  familiar  curl  curl  E_  -  <./|ikE    0,  using  a  steady  state  time  factor 
So] "'  '  "i'      to  this  using  the  standard  technique  would  next  be  related  to 
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scalar  Kelmholtz  solutions   However,  the  vectors  M  =  curl  A\y  and  N  =  1/k 
curl  curl  A$  where  \y  and'$  are  scalar  solutions  and  A  is  either  a  constant 
vector  or  a  unitary  vector,  are  no  longer  solutions  of  the  vector  equation, 
The  field  can  be  found  in  theory  in  terms  of  two  scalars  but  the  relations 
appear  unmanageable  for  this  case^^.   Direct  numerical  methods  are  a 
possibility  but  are  quite  difficult.   Second  is  the  problem  of  the  scalar- 
equation.   In  cylindrical  coordinates  after  the  time  and  z  variations 
have  been  cast  out,   the  scalar  wave  equation  becomes 

r  3,  (r  9wJ  +  a!|  +  (k2  -  32)  r2\|/  -   0,     k2  =  w'n  e,     (6  1) 

3r    3r    9cp 

3  is  the  z  coordinate  eigenvalue   The  crux  of  the  problem  is  to  separate 
this  two  dimensional  equation  so  that  a  solution  can  be  found.   It  may 
well  be  possible  to  do  this,  but  such  an  algorithm  has  not  yet  been  de- 
vised. 
6  2  Sectoral  Dielectric 

The  problem  of  a  continuous  dielectric  constant  variation  is  not 
exactly  soluble,  a  cogent  reason  for  considering  a  discreetly  inhomogen- 
eous  configuration.   Such  might  consist  of  sector-shaped  slabs  of  dielec 
trie  filling  the  radial  line,  each  sector  possessing  a  different  dielec- 
tric constant.   (See  Fig.  23.)   In  an  actual  case,  a  continuous  dielec- 
tric distribution  would  probably  be  approximated  by  a  step  function  of 
dielectric  sectors.   The  sector  problem  is  thus  a  very  practical  problem. 
Several  similar  problems  have  been  solved,  see  Marcuvitz  ^^'  and 
Klopfenstein^'*' .   The  method  of  attack  is  to  express  the  general  field 
in  one  region  as  a  sum  of  the  normal  modes  particular  to  that  region, 
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Figure  23   Sector  Dielectric  Loading 


and  to  do  this  also  in  a  contiguous  region   Then  the  tangential  E  and 
H  fields  are  equated  at  the  interface.   This  is  repeated  for  each 
discontinuity  and  the  resulting  equations  solved  for  the  coefficients. 
In  the  case  of  a  radial  or  coaxial  line  with  concentric  dielectrics, 
the  procedure  is  straightforward  but  lengthy,  although  the  solution 
of  a  transcendental  equation  may  be  required   An  additional  degree 
of  exacerbation  is  introduced  when  the  boundaries  are  radial  rather 
than  concentric   A  set  of  simultaneous  transcendental  equations 
arises  which  seem  to  defy  the  analyst   The  difference  between  the 
two  cases  arises  in  the  following  way,  where  a  radial  guide  with  no 
z  variation  is  used  for  an  example   The  two  eigenfunctions  are 
Bessel  functions  and  cosines  (also  sines  in  general)    For  the  con 
centric  boundary  when  the  tangential  fields  are  equated,  the  cosine 
factor  is  the  same,  term  by  term,  on  both  sides,  and  the  equality  must 
bold  for  all  'p   Thr>  Bessel  function  parts  are  equal,  term  by  term 
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Unconnected  equations  result  which  are  readily  solved   With  a  radial 
boundary,  on  the  other  hand,  the  cosine  factors  are  fixed  whereas  the 
Eessel  function  factors  differ  in  eigenvalue  and  cannot  be  individually 
equal;  perforce,  one  must  consider  the  series  in  toto   Here  then  is 
the  set  of  simultaneous  transcendental  equations 
6  3  Approximate  Solutions 

Since  the  problem  cannot  at  present  be  exactly  solved,,  approxi 
mate  methods  are  in  order.   One  such  method,  valid  for  coaxial  or 
radial  lines  whose  ratio  of  outer  to  inner  diameter  is  nearly  unity, 
has  been  obtained  by  considering  the  transmission  region  to  be  a 
rectangular  waveguide  wrapped  around  laterally ^5)    Another  attack 
is  that  of  Schelkunof f ^o;  wnich  replaces  the  Maxwell  equations  by  a 
set  of  ordinary  differential  equations.   If  the  parameters  of  the 
medium  can  be  linearly  approximated,  then  the  set  of  differential 
equations  reduces  to  a  set  of  linear  algebraic  equations,  a  situa- 
tion not  nearly  so  formidible„   It  would  probably  be  possible  to 
attain  a  good  answer  through  use  of  variational  methods  or  Rumsey's 
reaction  concept.   The  author  leaves  this  as  a  worthy  problem, 

A  crude  approximation  which  gives  usable  answers  assumes  the 
field  in  each  sector  to  be  that  which  would  exist  in  a  radial  line 
completely  filled  with  the  dielectric  of  that  sector.   This  is 
tantamount  to  neglecting  the  reflected  waves  at  the  discontinuity 
surfaces,  and  might  be  termed  an  optics  or  short  wavelength  approxi- 
mation  The  radial  guide  with  two  180°  sectors  of  dielectric  will 


exemplify  this  approach   Refer  to  Fig.  24. 
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Figure  2H       Two  Sector  Guide 

The   homogeneous    field    in   region   I   would   be,    from  Eq.    5-2 


Jzl 


ZTtC    Zjqq 


(6-2) 


with  the  current  I  flowing  in  the  feed  post  of  radius  c.   For  small  c 
(with  respect  to  A)  the  behavior  of  Z0o  will  be  inspected.   From 
Chapter  4, 


Z0  '     -    -  'TJToJi  -  Ni 


(6-3) 


1      2 

From  that  same  chapter,  W0    is  given  in  terms  of  S  ,  S    These  in  turn 
go  back  to  the  S(m,a),  T(m,a)  series  which  must  be  evaluated  for  small 
radius,  c 
With  kc  «1, 


h  (2)(kc)  ^   j  (2n)! 

2n  n!  (kc)n+1 


(rhn)^- 


j  (2n)! 


2n  (n-1)!  (kc)n+1 


(6-4) 


(6-5) 
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which  allows 


nc 


~-l  .  (6-6) 


(rhn)6 
When  these  are  inserted  in  Eq.  4  24,  the  series  reduce  to 

S(0,c)  =  -kc  So   ,   T(0,c)  =  0  . 

Thus  S  and  S  are,  respectively, 

S1  -  -kc  S0   ,   S2  =  0   .  (6-7) 

The  equations   for  U,    (Eqs.    4-9,    4-10)   now  give   the   value 

s'No  +  Ni 
-GTo1   =  ~ "Ob2  =  0  (6-8) 

S'Jq     +    Ji 

Using  Eq.  6-3  to  get  Z0o  in  terms  of  Bessel  functions  and  S's: 

S1  (J.No  -  JoNj 
Z-oo  ~  — (o-y) 

S'Jo  +  Ja 
For  small  argument, 

Jo  ^  1,   J!  ^Y2kc,      Nx  ~  ~  2 


jkc 


and 


Zq< 


S0  (4  +  TtkV  N0) 


n  kc  (2  S0  -  1) 

2   2 

For  wires  of  small  but  finite  size,  K   k  c  N0  <<;  4  and 

Zoo  -   4  S0Akc(4S0-l)   . 
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The  constant  part  may  be  called  zeta, 


4S0  C 

C  = ,  Z00  = •         (6-10) 

kc(4S0  -  1)  rckc 


The  electric  field  formula  becomes 

jrukjl  Zo(kjr) 


Ez2  "  ' 

2£  (6  11) 

Since  nk  =  wu.,  the  e  dependence  enters  only  in  the  Z0  argument.   The 
two  fields  are  then  given  by 

E     jri  k  I  ZoCkjr)     E  =  JT)  k  I  Z0(k2r) 

^zl  = "-"zl • 

2  t  2   C 

These  step  function  components  can  be  expanded  in  a  cosine  Fourier 

series  of  the  form 

oo 

ez  =  Y_  Zm(kr)  ^  °os  m  (3~i8) 

m=0 

with  the  coefficients 

A,  Zm  -  ^  M   Sin(^}   [Z0(kia)^(l-  [-1]-)  Z0(k2a)].   (6-12) 
mem£ 

6  ^  A  Twenty  Sector  Example. 

To  offer  an  example  that  could  yield  practical  answers,  one  can 
choose  any  large  number  of  sectors  in  which  to  divide  the  circle,   A 
convenient  number  is  twenty:  since  the  loading  was  to  be  symmetrical, 
thil   ivs  ten  variables.  (See  Fig  25.)   The  plane  of  symmetry  is 
represented  by  <p  -  0   For  each  sector,  the  field  is  given  (approxi- 
nately)  by  Eq   6  11    Again  the  step  function  may  be  expanded  in  a 
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certain    type   of  Fourier   serie.;: 

2jr\kKl 


JL 
10 


AmZm 


JL 

5 


*£ 


[  Z0(k!r)   J    cos  mcpdcp   +  Z0(k2r)   /  cos   mqpdqp 


JL 

10 


+  ...+   Z0(k10r)    i    cos   mcpdcp]     (6    13) 


_i2L 

10 


Figure  25   Radial  Line  With  Twenty  Sectors 


After  performing  the  indicated  integrations  the  result  is 


AmZm=-^ {ZoCkir)  [sin -sin  0]  + 

m   m     Cmem  10 


9  mil 


+  Z0(k10r)  [sin sin  mrt]  } 


10 


(6-14) 


Each  AmZm  is  thus  a  linear  combination  of  functions  of  the  ten  dielectric 

constants.   It  is  clear  that  if  the  ten  dielectric  constants  were 

given,  the  first  ten  Aj^Z  '  s  could  be  uniquely  found.   The  pattern  is 

approximated  by  the  terms  for  m  from  zero  to  ten. 

The  next  problem  is  to  find  the  radiated  or  far  field  corresponding 

to  these  coefficients.   The  azimuthal  field  is,  from  Chapter  2::;' 
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E^~nl   (rhn}'    VPA^O)    ^cos   m<p 
kr      "' 


+  H  mhn  Unm  Pnm(0)   Dnm   cos   mq> 


n  m 


(6-15) 


where  the  Anm,  Dnm  are  related  to  the  P^   by  Eqs .  2-14,  3-20,  amd  3-33 
Insertion  of  these  into  Eq.  6-15  yields,  after  rearrangement, 


(Unm)   cos  mcp 

2  2  - 

n  m   n(n  +  1) 


m2hn  [Pnm(0)j 


na 


y„ 


f  ka  (rhn)'    , 

[Pnm(0)] 

kr  (rhn)' 


(6-16) 


This  summation  may,  like  the  summation  of  Chapter  4,  be  rewritten  with 
the  m  sum  outermost.   Then  the  n  sum,  a  function  of  m  and  a  only,  is 


called  W   for  convenience,, 

»    (Unm)2  fka  (rhn)' 

w  -  y  — ^— — —  tpnm(o)]2- 

m   Z^*  n(n  +  1)1  kr  (rhn)  'a 


n-m 


m2hn[Pnm(0)]2 


na 


.     V 


\   Zm  Wm   cos   m® 


(6-17) 


(6-18) 


m=0 


The  radiation  field  (that  field  at  great  distance)  is  the  only  field  of 
interest  here  soWm  can  be  simplified  by  using  large  argument  approxi- 
mations for  the  Bessel  functions.   (See  Eq.  5  18).   In  the  form  given 
below,  the  constant  factor  kr  eJ    has  been  deleted. 


*m    '   l_ 


(Unm)2    (-j)n   fka[Pnm(0)]2  m2[Pnm(0)] 

+  J" 
n~~m*    r'(n  +  1) 


(6-19) 


(rhn); 


na 


Wm  h.i  ,  been  computed  using  techniques  already  described,  lor  several  ka 
and  is  shown  in  the  set  of  graphs,  Figs.  26  through  32..   For  a  given 
antenna  size,  ka,  and  for  a  given  set  of  dielectric  sectors,  the  AmZm 
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may  be  found  from  Eq  6-14,  the  Wm  from  the  curves;  these  values  are 
inserted  into  Eq.  6- 18  to  get  the  radiated  field  pattern.   This  completes 
the  approximate  analysis. 
6.5  Synthesis  Methods, 

If  the  phase  and  magnitude  of  the  desired  pattern  is  specified,  the 
problem  is  comparatively  simple;  the  pattern  is  expanded  in  a  Fourier- 
Legendre  series  and  like  coefficients  are  matched.   Usually,  however, 

only  the  magnitude  of  the  desired  pattern  is  known  and  there  is  no  way 
of  deciding  what  the  most  suitable  phase  variation  is.   This  synthesis 
problem  has  no  unique  solution,,  i.e.,  more  than  one  set  of  coefficients 
will  produce  the  same  amplitude  distribution. 

Several  methods  of  attack  have  been  used  for  linear  arrays,   Eaton, 
Eyges,  and  Macf  arlane^ '  '    work  with  the  absolute  value  of  the  space 
factor*,  expanding  the  square  of  the  absolute  value  in  a  Fourier  series 
The  solution  does  not  seem  possible  for  other  than  polynomial  formulations 
Another  method,  that  of  Taylor  and  Whinnery ' *"'    is  similar  to,  but  dif- 
ferent from,  that  of  Eaton  et.  al    It  involves  adjusting  the  space 
factor  to  make  it  an  analytic  function,  then,  since  the  magnitude  and 
phase  of  such  functions  are  intimately  related,  it  is  theoretically  pos- 
sible to  determine  the  phase  associated  with  the  magnitude  function. 
Implementation  of  this  for  linear  arrays  is  quite  difficult  and  appears 
not  feasible  for  an  equation  of  the  form  of  Eq  6  18   The  remaining 
■tteck  appears  to  be  that  of  assuming  a  constant  phase  pattern 


•The  space  factor  is  essentially  the  angular  variation  of  the  array 
pattern  i n  a  series  form 
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where  lobes  of  the  pattern  are  either  in  phase  or  out  of  phase  by  180 

degrees   This  is  accomplished  by  expanding  the  desired  pattern  U(cp) 
in  a  truncated  Fourier  series,  then  choosing  the  coefficients  of  Eq. 
6-18  to  make  each  term  of  that  equation  real  and  equal  to  the  corres- 
ponding term  of  the  series  for  U(cp) 

Let  U(qp)  be  a  symmetrical  pattern  with  the  axes  chosen  so  that  t lie 
Fourier  expansion  is  in  cosines  alone, 

00 

' —  e     2K 

U(qp)    =     '       bm  cos  mqp        ,      bm   =  _JH.  J     U(<p)    cos  mcpdcp,  (6-20) 

fi~0  2n     ° 

Under   the   conditions   of   the   above   hypothesis, 

K-  Re   (\  Z,  Wm)      ,      0  -   Im   (A^WJ      .  (6-21) 

After  inverting  these  for  A^^,  one  has 

bm  R«  <WJ                 "bm  I"  <V        ,     , 
Re(AA)— J_,   fa(W.__ .     (6-22, 

It  is  desired  to  equate  these  to  the  Aj^Z^j  in  terms  of  the  sector  parame-ti 
ters   (See  Eq.  6-14).   From  an  earlier  portion  of  this  paper,  ZQ(ka)  = 
1ST1  Jo  +  N0)  +  J^Jq.   This  allows  Eq..  6-14  to  be  split  up  into  real  and 
imaginary  parts,  and  equated  to  the  forms  above,  with  the  result 


b     Re    (WJ  t  -1  mn 

_= {[•5r2j0(k1a)]     [sin sin  0]    +••••} 


x 


AM'  2nknl     mem  10 


bm  Im   (WJ  C  -1         M  tan 

-    x  = {[TJj0(k1a)    +  No(kia)]     [sin -sin  0]    +    •••} 


|W    r  2nkrtl     mem  10 

m  (6-23) 
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By  inspection  it  is  seen  that  the  J0  terms  of  the  upper  and  lower  equations 
are  proportional;  these  may  be  eliminated  between  the  two  to  produce  a 

single  set  of  equations  relating  the  sector  parameters  to  the  desired  field. 

Only  the  relative  magnitude  is  of  importance  so  the  constants  on  the  left 

hand  side  have  been  replaced  by  a  single  constant 


Const,    x   bm    [i/lm  Wm  -  TD^Re  Wj 


IW    I2 
1    m > 


-1 


mc. 


ran                                                                   9mn 
(Not^a)  [sin -sin  0]    +    •  *  •+  N0(kl0a)     [sin sin  mrc]  } 


10 


10 


(6-24) 
It  will  be  necessary  to  adjust  the  constant  to  scale  the  solutions  into 
the  range  of  the  No(kjLa)  function  which,  for  all  practical  cases,  will  be 
less  than  one  in  absolute  value.   The  £•  or  |i^  are  found  from  the  k^  need- 
ed to  fit  the  NQ  to  the  values  found  from  the  Eq.  6-24   this  set  of  equa- 
tions must  be  solved  simultaneously  for  the  ten  N0   With  high  speed  com- 
putational aids  available,  this  is  exceedingly  simple;  without  such  help 
the  synthesis  technique  may  not  be  possible   The  bracketed  sine  coeffi- 
cients of  Kq   6-24  are  listed  in  TabJe  8 

Some  latitude  exists  in  the  choice  of  antenna  radius  a   As  in  all 
types  of  arrays,  a  given  pattern  will  require  more  complex  excitation 
(widely  varying  dielectric  constant)  as  a  decreases,  resulting  eventually 
in  a  supergain  phenomenon  where  rapid  large  changes  in  e  versus  qp  take 
placi    A  ten  sector  configuration  would  not  approximate  such  a  situation 
Probably  a  good  starting  value  would  be  the  radius  of  a  circular 
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m  =  3 


.15451 

. 14695 

.13483 

.13938 

09082 

.02367 

.11061 

00000 

-.10701 

.07102 

-.09082 

-. 14947 

.02447 

-. 14695 

-06870 

.02447 

-.14695 

.06870 

.07102 

-.09082 

14947 

.11061 

.00000 

.10701 

.13938 

.09082 

-.02367 

15451 

.14695 

-.13483 

11888 

. 10000 

.07925 

04541 

-. 10000 

-.12824 

14695 

-.10000 

.00000 

04541 

. 10000 

. 12824 

11888 

. 10000 

-.07925 

11888 

-.10000 

-.07925 

04541 

-. 10000 

. 12824 

14695 

. 10000 

. 00000 

04541 

. 10000 

-.12824 

11888 

-.10000 

.07925 

m  -  7 


m  - 


m  -  9 


.05779 

.03674 

.01717 

.12572 

-.09618 

-.04982 

.09000 

.11888 

.07760 

.01991 

-.09618 

-.09778 

.11341 

.03673 

.10839 

.11341 

.03673 

-.10839 

.01991 

-.09618 

09778 

.09000 

.11888 

-.07760 

.12572 

-.09618 

.04982 

.05779 

.03673 

-.01716 

Table  8. 
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array  which  produces  a  comparable  pattern^4' 


Two  examples  will  be  worked  out.   The  first  is  a  cardioid,  the  sec 
ond  an  optimum  fan  heam 
6  6  The  Cardioid. 

The  cardioid  pattern,  |E|  -  1  +  cos  cp,  offers  an  elementary  example. 
It  is  already  written  as  a  Fourier  series  of  two  terms,  so  immediately 


b0  =  1,   bi  -  1,   b   -  0   for  m  >  1 


Equation  6  24  becomes 


-V    ImV.'0   -VfieWo  *  ic 
=  -{N0(kia)  —  +   ••-    +  N0(k10a) — } 

|W0'2  10  10 


■ctq    Im  Ws   -  T«J"o  Re  W  i  n 
=    -%{N0(kaa)     [sin  —  -sin   0]    + 

Iw-J2  10 


9TC 

+   IM0(k10a)    [sin -sin  Tt]  } 

10 


27t  18Tt 

0   a    {No(k.a)    [sin sin  0]    +••«    +  N0(k10a)    [sin sin  2k] } 

10  10 

(6-25) 
The  desired  pattern  has  been  analyzed  into  Fourier  cosine  coefficients 
and  these  have  been  inserted  into  Eq  6  24   Next  a  normalized  radius, 

selected  such  that  the  decrease  of  Wm  with  increasing  m  is  commen 
with  the  corresponding  decrease  of  the  bm    In  this  example, 


85 
ka  -  2  83  was  chosen.   Values  of  VV   are  then  put  into  Eq   6  25  and  the 
set  solved  with  the  answers  of 


N0(kia)  -  const  x  -  1106 

N0(k2a)  =  const  x  -  0975 

N0(k3a)-  const.  x  -.0727 

N0(k4a)  -  const,  x  -.0386 

N0(k5a)  -  const,  x  0014 

N0(k6a)  =  const  x  0436 

N0(k7a)  =  const,,  x  0837 

N0(kea)  =  const  x  H78 

N0(k9a)  =  const ,  x  1426 

N0ik10a)  -  const,  x  1556              (6-26) 


A  suitable  scale  factor  or  constant  is  3,81  The  values  calculated  for 
6^  are  listed  in  Table  9-  It  would  be  possible  to  produce  sectors  with 
the  proper  dielectric  constant  by  loading  paraffin  with  titanium  dioxide 


Sector  Solutions  k^a  NQ(k^a)  /e^  e^ 

i  to 
6  25 

1  -.1106  4  89  -290  1  73  3  00 

2  -.0975  4,73  -.256  1.67  2, 80 

3  -.0727  4  49  -191  1.58  2.50 

4  -  0386  4  22  -.101  1  49  2  22 

5  0014  3  95  0037  1  39  1  94 

6  0436  3  68  115  130  1  69 

7  0837  3  42  220  121  1  47 

8  1178  3  19  309  1  13  1  28 

9  .1426  3  01  374  1  06  1  12 
10  .1556  2  90  408  1  02  1,04 

Ta^le  9   Dielectric  Sectors  for  Cardioid 


powder;  or  by  perforating  a  suitable  material  with  an  array  of  small  holes 
to  adjust  the  dielectric  constant   However,  satisfactory  results  can  be 
obtained  from  a  simpler  arrangement    Measurements  were  taken  (see  Chapter 
7  for  a  description  of  the  experimental  apparatus)  using  a  180°  sector  of 
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polystyrene,  dielectric  constant  of  2  5   This  approximates  five  of  the 
sectors  of  Table  9  by  2  5  and  the  remaining  five  sectors  by  10   The 
average  dielectric  constant  was  then  1  75  which  made  the  average  ka  equal 


to  i  1  75  ka  or  3  75   This  corresponded  to  a  model  frequency  of  2340  mc 

In  Figure  33  is  depicted  the  theoretical  curve  together  with  the  measured 

curve   Excellent  agreement  is  noted  except  in  the  rear  direction;  prob- 
ably a  more  accurate  mockup  of  the  sector  dielectric  constants  would  improve 

this  condition   The  6  pattern,  for  qp  -  0  is  shown  in  Figure  34 
6  7  Fan  Beams. 

For  an  example  of  shaped  beams,  an  optimum  fan  beam  with  equal  level 
side  lobes  will  be  worked  out  This  pattern  is  a  Tchebyscheff  polynomial 
in  the  azimuthal  plane  and  a  fan  shape  in  elevation  Much  space  has  been 
devoted  in  the  literature  to  array  synthesis  and  to  the  use  of  Tchebyscheff 
polynomials  to  optimize  certain  patterns  References  4  and  38  contain  ex- 
tensive and  up  to-date  bibliographies  Tchebyscheff  polynomials  are  de 
fined  by(39) 

Tj\(z)  -  cos  (N  arc  cos  -  )      ,       \  z  |  <_  1 

TN(z)  --  cosh  (N  arc  cosh  z)     |z|  >  1      (6-27) 

where  it  is  convenient  to  let  z   c  cos  cp  +  d   The  amplitude  of  the  func- 
t  ion  varies  from  H  to   1  and  back  to  H  as  cp  goes  from  0  to  2k  Oscillation 

•  between  the  range  +  1  which  gives  the  function  a  "main  beam" 
plus  ,i  ouaber  of  equal  level  "side  lobes"   For  a  given  side  lobe  level  1/R 
"1  >'  the  ratio  of  main  beam  to  side  lobes)  the  range  of  argument  z0  is  found 


^ 


V 


/ 


/ 


magnitude  of  Eq  is  plotted  versus 
solid  line:   theoretical  pattern 
dashed  line:   measured  pattern 


Figure  33.  Cardioid  Patterns  of  Spherical  Antenna 


agnitude  of   £q   is    plotted  versus  9 


Figure  34.     Elevation  Pattern  of  Spherical   Antenna 
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from  z0  ~  cosh  ( 1/N  arc  cosh  R);  c  and  d  are  found  from 

c  -  %   Uo  +  1)   ,  d  -  %  (z0   -  1)   .        (6  28) 

The  Tchebyscheff  polynomials  can  be  expanded  in  an  exact  finite  Fourier 
series^''  of  the  form 

N 

TN  (cp)  =  Y   CNX  cos  icp  .  (6-29) 

Duliamel  gives  the  coefficients  for  N  up  to  eight. 
For  the  example,  take  N  =  5.   The  C's  are 

C05  ■-   5d  -  20d3  +  16d5  -  30c2d  +  8c2d3  +  30c^d 
Ci5  =  5c  -  60cd2  +  80cd4  -  15c3  +  120c3d2  +  10c5 

C26  =  -30c2d  +  80c2d3  +  40c4d 

C35  -  -5c3  +  40c3d2  +  5c5 

C45  =-    10c4d 

C56=  c5  (6-30) 

Let  R  -  5  which  corresponds  to  -14  db  side  lobes.   Then 

z0  =  1-10695898 
c  =  1.05347949 
d  -  0.05347949 

and  the  values  calculated  for  the  series  coefficients,  calling  them  bm 
now,  are  the  following- 
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b0  =  0.46124418  bt  =  0.92687461 
b2  =  0.86782692  b3  =  0.77574232 
b4  =  0.65870596        b5  ■  1.29756889-        (6-31) 

The  pattern  calculated  from  these  b's  is  that  of  Figure  35  where  the  main 
beam  to  side  lobe  ratio  of  five  is  seen.   The  b  values  were  inserted  into 
Eq.  6-24  along  with  the  functions  for  ka  =  5.80  which  was  taken  for  a  typical 
example.   The  answers  and  the  calculations  for  £j  are  given  in  Table  10. 
Again  rather  than  attempting  to  produce  the  exact  dielectric  constants, 
polystyrene  and  air  were  used.   It  was  thought  that  the  values  of  Table  10 
might  be  approximated  by  four  sectors  of  polystyrene,  five  of  air,  and  one 
of  polystyrene,  giving  again  an  average  dielectric  constant  of  1.75  with 
the  corresponding  average  ka  of  7.66  and  model  frequency  of  4800  mc.   Since 
polystyrene  has  a  dielectric  constant  considerably  higher  than  that  needed 
here,  smaller  sectors  (in  angle)  are  indicated.   Measured  patterns,  using 
seven  sectors  of  polystyrene  in  the  cp  =  0°  direction  and  two  sectors  of 
polystyrene  in  the  cp  :  180°  out  of  a  total  of  20  sectors,  are  shown  in 
Figs  36  and  37.   The  beam  has  a  ratio  of  4  8  which  is  very  satisfactory. 
Additional  patterns  for  different  elevation  angles  are  also  shown.. 

It  is  interesting  to  note  that  the  theory  demands  dielectric  loading 
in  both  the  main  beam  direction  and  in  the  rear  direction  and  that  this 
requirement  is  borne  out  by  the  measurements., 
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Sector 

Solutions 

kia 

N0(kia) 

i 

to 
Eq.  6-24 

1 

1609 

7  86 

199 

1.36 

1.85 

2 

.0972 

7,35 

079 

1  27 

1  62 

3 

0594 

7.25 

048 

1.25 

1  56 

4 

0883 

7.33 

071 

1.26 

1.59 

5 

0972 

7.35 

.079 

1.27 

1  62 

6 

-.0144 

6,69 

-018 

1.15 

1.33 

7 

-.1927 

6.25 

-.238 

1.08 

1.17 

8 

-.2566 

5.  80 

-.318 

1  00 

100 

9 

-. 1270 

6  39 

-  103 

1  10 

1,21 

10 

.0428 

7.27 

.053 

1.25 

1.56 

Table  10   Dielectric  Sectors  for  Fan  Beam 


Tchebyscheff   T5   polynomial    plotted 
versus   cp 


Figure  35.  Tchebysct.eff  Fan  Beam 


curve  is  measured  value  of  magnitude  of  En 
versus  9;  <p  =  0. 


Figure  37   Elevation  Pattern  of  Fan  Beam  Antenna 
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9=75° 


G=60° 


6  =  45° 


6=30° 


cp  patterns        6  polarization 

Figure  38.   Fan  Beam  for  Different  Elevation  Angles 


7.   THE  EXPERIMENTAL  FACILITY 

7  I  The  Swing  Frame  Ground  Screen 

Considerable  effort  has  been  expended  by  the  author  and  Roger  R. 
Trapp  in  the  design  and  construction  of  equipment  to  measure  principal 
plane  patterns  of  monopole  antennas  or  of  other  antennas  requiring  a 
ground  plane,  e.g.,  slot  antennas.   The  setup  consists  of  a  swing  frame 
ground  screen  with  rotating  aperture,  together  with  the  necessary  elec- 
tronic equipment.   Frequency  coverage  is  1  to  18  kmc;  at  1  kmc  the  screen 
is  12  wavelengths  on  a  side. 

The  screen  is  12  feet  square  and  is  supported  by  a  structural  steel 
and  wood  framework  five  feet  high.   This  height  allows  the  operators  to 
see  over  the  screen,  a  factor  that  is  more  important  than  large  height 
under  the  screen.   The  location  is  the  roof  of  the  Electrical  Engineering 
Building,  adjacent  to  an  equipment  shack.   "Boiler  Plate"  construction  was 
used  for  the  surface  of  the  plane;  4"  by  6"  sheets  of  dural  %  inch  thick 
were  lapped-joined  to  provide  a  smooth  conducting  surface.   (See  Fig.  39.) 
Joints  were  first  milled  and  then  hand  scraped,  with  a  screw  spacing  of 
one  inch  used  to  fasten  the  joints  together.   The  entire  top  surface  is 
rigidly  bolted  to  the  framework  only  near  the  center  to  allow  expansion 
due  to  the  wide  variation  of  temperature  in  this  climate  (about  50°  C.  )  , 
In  the  center,  segmented  phosphor  bronze  sliding  fingers  are  mortised 
into  a  1/32  inch  deep  counterbore  on  a  26  inch  circle  in  the  aluminum 
plate   These  stationary  fingers  slide  on  an  outside  counterbore  in  a 
rotating  brass  ring;  this  ring  contains  a  circular  aperture  plate  of  21 
iriffi  diametei    Support,  for  the  ring  is  provided  by  an  aluminum  anulus  of 
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24  inch  0   D, ,  two  inches  thick  and  two  inches  deep,  which  is  secured 
to  the  underside  of  the  brass  ring.   The  aluminum  anulus  is  angle 
ironed  to  a  B-29  gun  turret  ring  and  bearing  assembly.   The  latter  has 
two  concentric  formed  rings,  about  a  yard  in  diameter,  with  one  ring 
supported  from  the  other  by  a  number  of  radial  and  thrust  ball  bearing 
races   Figure  40  shows  the  underside  of  the  screen,  looking  up  at  the 
ring  assembly   This  type  of  construction  allows  the  access  space  be- 
low the  aperture  plate  to  be  completely  clear   Drive  is  provided  by  a 
General  Radio  1/15  hp  direct  current  motor  through  a  Boston  gear  re 
duction  box,  using  the  large  ring  gear  of  the  gun  turret  assembly  for 
traction   Position  data  are  provided  by  a  1:1  geared  Selsyn. 

Another  aperture  plate  contains  a  square  recess  to  accommodate 
aperture  plates  14  inches  square,   Both  plates  fasten  by  flat  head 
screws  about  1  1/4  inch  apart.   The  top  surface  of  the  screen  and  aper 
ture  is  flat  and  smooth 

A  boom  or  swing  frame  carries  the  receiving  horn.   As  may  be  noted 
from  the  photograph,  the  boom  is  pivoted  at  two  gear  boxes  located  on 
opposite  corners  of  the  screen.   Vertical  members  are  2"  by  4"  lumber 
while  the  horizontal  crosspiece  is  a  composite  hollow  four  inch   square, 
glued  and  dowelled   A  universal  mounting  bracket  allows  quick  inter 
changeability  of  antennas   The  gear  boxes  each  contain  a  large  worm 

•r  to  which  the  2  by  4  is  bolted,  and  a  worm  drive;  one  box  also 
contains  a  Selsyn  takeoff   These  gear  boxes  are  connected  by  a  shaft 
train  which  La  offset  to  avoid  fouling  the  aperture  access,  and  which  is 
'irivr-n  near  its  center  Ijy  another  DC  motor  through  a  gear  reduction  box. 
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Ball  bearing  pillow  blocks  were  used  throughout  to  reduce  friction  To 
insure  reliable  operation  in  damp  and  cold  weather,  the  motors  and 
Selsyns  were  encased  in  sheet  aluminum  boxes,  each  heated  by  a  light 
bulb,   Adjustable  limit  switches  prevent  the  boom  and  antenna  from 
hitting  the  ground  plane   The  crosspiece  is  adjustable  in  distance 
from  the  aperture  in  six  inch  steps  with  a  minimum  of  four  feet  and  a 
maximum  of  seven  feet   The  control  equipment  will  next  be  described. 

7  2  Control  and  Electronic  Equipment 

Power  for  the  drive  motors  is  supplied  by  a  five  ampere   Accurate 
Engineering  Co,  dry  rectifier  for  the  armatures  and  a  one  ampere  supply 
of  the  same  type  for  the  fields   These  supplies  are  also  used  by 
another  pattern  range   A  main  control  panel  in  the  equipment  shack 
(see  Fig   41)  contains  speed  adjustments  (variable  resistance  in  arma- 
ture circuit)  and  control  switches.   In  this  photograph,  the  racks  on 
the  left  hand  side  belong  to  another  pattern  range   Motors  may  be 
started  and  reversed  from  three  locations  main  control  panel,  recorder 
console,  and  a  switchbox  on  the  ground  screen,  with  relays  providing 
the  necessary  interlocking   Selsyn  information  from  the  ring  or  boom 
may  be  switched  to  a  dial  indicator,  or  to  the  polar  recorder   A 
Styroflex  (low  loss)  cable  and  several  other  coaxial  cables  connect 
the  screen  and  the  equipment  racks   These  racks  contain  several 

illators  of  the  cavity  type  covering  the  range  300  to  3000  mc  as 
well    ;is  Klystron  power  supplies  used  to  pipe  voltages  over  a  special 
cable  to  Klystrons  that  might  be  located  at  the  antenna  terminals. 
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To  date  it  has  been  possible  to  cover  the  entire  range  with  adequate  power 
using  the  Hewlett-Packard  line  of  standard  Klystron  signal  generators.   A 
versatile  switching  arrangement  allows  any  oscillator  to  be  connected  to  a 
wavemeter  or  to  the  aperture. 

The  received  signal  is  intercepted  by  one  of  a  set  of  horns  and  corner 
reflectors  that  cover  the  range  1  to  18  kmc .   A  crystal  or  bolometer  de 
tector  may  be  used,  feeding  the  1000  cycle  signal  to  a  PRD  bolometer  ampli- 
fier via  a  coaxial  switch  panel .   A  signal  is  taken  from  this  amplifier  to 
operate  an  external  portable  meter,  and  also  the  recorder.   This  latter 
signal  is  rectified  by  a  synchronously-driven  chopper  and  the  resultant 
DC  is  fed  to  the  recorder  circuits.   Both  linear  and  square  root  response 
are  available  in  this  recorder  (Leeds  and  Northrup),  an  unusual  feature 
of  which  is  that  the  servo  amplifier  gain  is  adjusted  (on  square  root)  by 
a  potentiometer  ganged  to  the  pen  shaft  so  as  to  maintain  the  servo  loop 
gain  at  the  low  end  of  the  scale   This  reduces  the  sluggishness  of  the 
pen  near  zero. 

To  obtain  patterns  from  an  antenna  mounted  on  either  the  square  or 
round  aperture  plates,  the  procedure  is  simple.   The  plate  is  bolted  in, 
then  switches  are  set  to  connect  a  suitable  oscillator.   Next  the  proper 
horn  and  detector  is  mounted  on  the  boom.   Maximum  signal  is  then  obtained 
tuning  the  matching  stubs,  using  the  external  meter.   The  aperture  may  be 
set  to  a  desired  <p  angle  and  the  boom  rotated  through  9,  or  the  boom  may 

et  to  a  particular  0  and  a  cp  pattern  taken   The  polar  recorder  chart 
t;ihle  servo  system  follows  in  synchronism  of  course. 


8 .   CONCLUSIONS 

A  technique  for  finding  the  radiated  fields  of  a  non  uniformly 
excited  spherical  dipole  has  been  demonstrated.   By  assuming  a  constant 
distribution  of  tangential  electric  field  across  the  gap,  a  method  is 
established  for  deducing  the  cp  distribution  of  the  gap  field.   The  same 
type  of  attack  will  not  suffice  for  spheroidal  or  ellipsoidal  dipoles 
as  the  simple  Hertzian  vector  solutions  do  not  exist  for  these  cases. 

An  offset  feed  can  produce  a  wide  variety  of  patterns  for  dif- 
ferent feed  distances  and  antenna  diameters.    Since  the  spherical 
antenna  depends  almost  solely  on  lengths  and  very  little  on  angles,  it 
is  frequency  sensitive,  hence  narrow  band.   In  the  plane  of  the  horizon, 
only  one  polarization  exists. 

The  synthesis  technique  developed  allows  azimuthal  patterns  (which 
are  admissible)  to  be  produced  by  sectoral  dielectric  loading.   In  the 
cardioid  and  fan  beam  examples  worked  out,  good  correlation  was  obtained 
with  measured  results,   A  phenomenon  similar  to  super-gain  occurs  in  that 
an  antenna  diameter  too  small  for  a  given  pattern  requires  large  ranges 
of  dielectric  constant,  and,  for  these  cases,  the  approximations  made  in 
the  solution  are  not  valid.   A  rotatable  beam  antenna  would  be  possible 
using  ferrite  material  and  an  electrically  rotated  magnetic  field  struc- 
ture inside  the  antenna. 
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APPENDIX    I 


CALCULATION   OF  THE   Sm   AND  Tm   SERIES 


The  Sm  series   was 


00 

S      =    Y      X2n    *    l)(n  -   ajj 

m       Z_*2n2(n   +   l)(n  +  m) ! 
n=m 


3Pr 


30 


3Pr 


39 


IL-A 
2 


(1-1) 


where  A  was   the   gap   angle.      For   argument   9   =    n/2,    the   Legendre    function 
derivative   has   a   simple    form: 


9pnm 


30 
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n+3rn  +  l 


1-1 -3° 5...     (n   +  m) 
1-2-4-6. ..     (n  -  m  - 


1) 


for    (n    +  m)    even 


for    (n    +  m)    odd 


(1-2) 


Thus    the   summation   is   over  odd  or   even   n   only,    depending  on   m.      The 
asymptotic   expansion   was   given    in    the   text    (Eq„    4-20);    this   was   obtained 
by   differentiating   the   given   series   with   respect    to    9   and  keeping  only 
terms   of   the   highest   order. 

nm+^  sin[(n   4   %)    Q  +  m*  -  tt] 
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4sin   G 
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The   approximation   with   respect   to  n    is   poor    for  n   as    large   as   10 

whereas    the   0  approximation   is   good    for   9  near  rc/2.      A   good   approximation 

for    both    is    gained    by   combining  Eq      1-2   and   Eq.    T-3    to   get 

p  m  (I-*) 

±I!_  -v   (    1)^1        M35 (n    4  m)  in    [(n   +  1/)Q  +  m  ..  aj 

12   4-6...    (n    -  m    -  1)  2        4 
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A  typical  case  would  have  a  gap  angle  of  1°  or  A  =  1°.   Then 


n+3m+l 
Pnm'(A)  -  (-1)   2      l|3-5    (n  +  m)    cQs  (n  +  %)Q 

1  2*4"6. . .  (n  -  m  -  1) 


(1-5) 


The  general  term  of  the  series  is  then 


(2n  +  l)(n  -  m)! 
2n2(n  +  l)(n  +  m) ! 


1-1'3-S...  (n  +  m) 


1-2-4-6.. .  (n  -  m  -  1) 


cos  (n  +  y2)°  (1-6) 


where  (n  +  m)  is  odd.   For  large  n  this  becomes  ^  cos  n.   If  a  zero  gap 

had  been  assumed,  the  cos  n  factor  would  be  absent  and  the  series  would 

diverge.   The  cos  nA  factor  produces  a  slow  sign  change  or  alternation; 

clearly,  the  larger  the  gap  the  faster  the  sign  change.   A  one  degree  gap 

corresponds  to  a  sign  change  every  180  terms. 

Since  (n  +  m)  is  odd>  two  cases  may  be  treated. 

Case    1 

n  odd,  m  even. 
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!  a  s  e  2 
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F„  =  (n  -  2)(n  t   D(n  +  3)    etr 
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The  general  scheme  for  performing  the  calculation  is  to:   1)  set  up  a 
table  of  cos  (n  +  lA)    in  the  memory,  2)  form  an  iterative  loop  for  com- 
puting 

'  1-1-3-5. ..  (n  -  2) 
1-2-4-6...  (n  -  1) 

using  previous  values.  3)  set  m,  4)  form  a  loop  for  finding  Fm  for  the 
given  m,  5)  change  m  and  repeat.   The  program  for  one  of  these  calcu- 
lations is  shown  in  the  next  several  pages  and  in  the  block  diagram  of 

Fig.  II. 

There  is  a  possibility  of  sign  change  every  90  terms;  the  signs 
alternate  in  pairs,  +,-,-, +,+,-,  etc.   The  terms  have  been  combined  in 
sets  of  two,  one  plus,  one  minus.   Each  set  then  contains  180  terms. 
The  sets  form  an  alternating  series  so  the  truncation  error  can  be  re- 
duced by  adding  only  half  of  the  last  set.   Using  the  ILLIAC,  five  sets 
(900  terms)  were  evaluated  giving  an  accuracy  of  four  significant  fig- 
ures with  the  Las!  figure  in  doubt.   Cumulative  error  is  believed  to  be 

than  the  truncation  error 
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The  code  prints  the  first  25  terms  (for  later  use)  and  the  partial 
sums  for  each  set,   Table  1-1  shows  the  results. 


m 


'in 


0 

1.6276 

0 

1 

1.0558 

1.1362 

2 

.8251 

.9991 

3 

.6868 

.9028 

4 

5896 

.8283 

5 

5150 

.7677 

6 

.4549 

.7163 

7 

.4046 

.6719 

8 

.3616 

.6326 

9 

.3241 

.5973 

10 

.2910 

.5653 

Table  1-1. 

Values  of 

Sm 

and  Tm 

Tm  Series 

The  T  series  was 
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Following  the  same  steps  used  before,  a  good  approximation  for  both  n  and 
6  is,  for  a  one  degree  gap 
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The   general    term   is,    where    (n    +   m)    must   be   even. 
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This  has  been  separated  into  two  cases  and  computed  in  the  same  fashion  as 
the  earlier  series.   Results  of  this  are  also  shown  in  Table  1-1. 


Program  for  Sm  Series  -  m  even. 
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0010K 
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term  code 
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L4 

L 

42 

31L 

2 

L5 

91F 

L0 

97F 

3 

32 

4L 

41 

92F 

_4J6 
L5 

28L 
91F 

5 

L4 

80F 

40 

2F 

6 

19 

10F 

L4 

91F 

7 

66 

2F 

7J 

95F 

8 

40 

IF 

50 

IF 

9 

7J 

95F 

40 

92F 

10 

41 

96F 

L5 

91F 

11 

40 

F 

LO 

80F 

12 

40 

IF 

L5 

2F 

13 

L4 

80F 

40 
14  L5 

3F 

96F 

LO 

97F 

15 

36 

28L 

L5 

96F 

16 

L4 

81F 

40 

96F 

17 

L5 

IF 

66 

F 

18 

7J 

92F 

40 

92F 

19 

L5 

3F 

10 

IF 

20 

66 

2F 

7J 

92F 

21 

00 

IF 

26 

22L 

22 

40 

92F 

L5 

F 

23 

LO 

81F 

40 

F 

24 

LO 

80F 

40 

IF 

25 

L5 

2F 

L4 

81F 

26 

40 

2F 

L4 

80F 

27 

40 

3F 

.26 

14L 

28 

L5 

91F 

L4 

80F 

29 

40 

IF 

LO 

80F 

30 

66 

IF 

7J 

95F 

31 

40 

95F 

22 

OF 

outer 

loop 

1 

30400K 

0 

92 

137F 

92 

512F 

1 

L5 

86F 

00 

9F 

2 

40 

86F 

L5 

82F 

3 

00 

9F 

40 

82F 

4 

L5 

98F 

00 

9F 

5 

40 

98F 

41 

97  F 

s 

'92 

259F 

92 

706F 

7 

92 

967F 

92 

706F 

8 

92 

194F 

92 

258F 

9 

92 

514F 

92 

194F 

10 

92 

706F 

92 

97  5F 

11 

92 

643F 

92 

967  F 

12 

92 

51 4F 

92 

706F 
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13  92  967F 
92  707F 

14  L5  97F 
10  29F 

15  JO  2F 
50  15L 

16  26  140F 
92  133F 

T792  513F 
50  17L 

18  26  10F 
92  145F 

19  92  513F 
26  20L 

20  L5  97F 
LO  98F 

21  32  23L 
L5  97F 

22  L4  81F 
40  97F 

23  26  6L 
OF  F 


tape  layout; 

1„  Decimal  Order  Input 

2.  0010K  mam  code 

3.  00100K  first  term  code 
4<,  00108K  general  term  code 

5.  00140K  P4  print  code 

6.  00170K  P2  print  code 

7.  00300K  Tl  sin-cos  code 

8.  00400K  outer  loop  code 

9.  24  999N  input  par.  tape 


parameter  tape: 

1.   0098K  OOmF  OOF 

2„   24  400N  start  program 


APPENDIX  II 

A  SHORT  SUMMARY  OF  THE  ILL  I  AC  ORDER  CODE 

The  machine  contains  an  accumulator,  A,  which  is  a  double 
shifting  register,  a  mul tipl ier -  quotient  register,  Q,  which  is  also 
a  shifting  register,  and  an  order  register.   Registers  A  and  Q  are 
coupled  together  so  that  a  number  may  be  shifted  any  number  of  places 

from  one  to  the  other.   In  addition,  one  starts  with  a  number  in  the 
accumulator  and  ends  the  same  way.   Multiplication  and  division  use 
both  registers   the  multiplier  is  placed  in  Q  and  the  product  (which 
has  79  digits)  appears  in  AQ  with  the  most  significant  portion  in 
A  the  dividend  is  placed  in  A,  the  quotient  appears  in  Q. 

The  input  to  the  machine  is  by  standard  five  hole  teletype  tape 
with  a  photoelectric  tape  reader.   The  sexadecimal  number  system 
is  used,  for  orders  and  numbers  on  the  tape    Output  is  by  the  same 
tape  via  a  high  speed  punch  or  by  a  cathode  ray  tube  console;  a 
special  program  is  provided  for  use  with  the  latter.   The  high 
speed  memory  contains  1024  locations  which  are  numbered  starting 
at  zero  with  each  spot  capable  of  storing  40  binary  digits    The 
registers  also  store  a  40  bit  number 

Orders  are  treated  in  the  machine  as  numbers,  hence,  arithmetic 
operations  may  be  performed  on  the  orders    This  fact  and  the  single 
;i'l<Jress  associated  with  each  order  put  the  ILLIAC  in  a  class  of 
machines  called  "Princeton  machines"  after  the  inventors,  von  Neumann 
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and  Goldstine   Each  order  consists  of  two  sexadecimal  instruction 
digits  and  an  address  number.   This  address  number  may  not  refer  to 
a  memory  location,  e.g.,  on  a  shift  order  the  address  gives  the  num- 
ber of  shifts   As  the  orders  occupy  less  than  20  bits  each,  two 
orders  are  packed  into  one  word  and  called  right  and  left  hand 
orders.   The  machine  executes  orders  in  sequence,  starting  at  a 
predetermined  point  in  the  memory:   LH  order  then  RH  order  of  first 
word,  LH  order  next  word,  etc.   The  sequencing  may  be  changed 
by  conditional  control  transfer  orders  and  stop  orders.   At  a  stop, 
two  possibilities  avail.   If  the  operator  uses  a  black  switch,  the 
machine  will  obey  the  stop  transfer.   If  the  white  switch  is  used, 
control  is  transferred  to  the  next  order.   Thus  the  course  of  a 
program  may  be  altered  by  the  operator. 

An  especially  useful  library  code,  the  Decimal  Order  Input, 
allows  the  orders  to  be  written  with  decimal  addresses  and  with  one 
of  six  terminating  symbols.   These  symbols  permit  maximum  flexi- 
bility in  the  use  of  subroutines  and  parameters.   They  are: 

K     when  OOpK  is  placed  on  the  tape,  orders  following  are 
stored  in  sequence  in  the  store  starting  at  location 
p,  and  p  is  added  to  all  relative  addresses. 

F     denotes  a  fixed  address. 

L     denotes  a  relative  address. 

J     for  input  of  occasional  decimal  fractions. 

N     transfers  control  to  program  at  location  q  when  used 
as  24  qN. 

S     for  preset  parameters,  S3  after  an  address  adds  the 
contents  of  location  3  to  the  address. 
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The    ILLIAC  coding   manual    provides    complete    and    detailed    infor 

(40) 
mation   on   the   machine   and    on    coding    methods    and    tricks 

The    following   is   an    abbreviated    order    list,     containing    the 

most    used    orders: 

shift   AQ   left    n    places,     l£n£63 

final    stop 

shift  AQ  right  n  places,  l£n£63 

make  A  ~  %  and  do  right  shift 

stop,  black  switch  transfers  control  to  RH  order  at  n 

transfer  control  to  RH  order  at  location  n 

same  as  20  except  LK  order 

same  as  22  except  LH  order 

if  A  >  0,  same  as  corresponding  2  orders 

if  A  ^  0,  go  on  to  next  order 

replace  number  in  memory  at  n  by  number  in  A 

replace  location  n  and  A  by  zero 

replace  address  of  RH  order  at  n  by  address  of  RII 
order  in  A 

same  as  42  except  both  LH  addresses 

replace  Q  by  number  in  n,  number  in  n  is  N(n) 

divide  A  +  2  39Q  by  N(n) 

put   N(n)    times    Q    +    2~39A    into    AQ 

multiply   N(n)     by    Q 

rounded    product    of    q    and    N(n) 


00 

n 

OF 

0 

10 

n 

19 

n 

20 

n 

22 

n 

24 

n 

26 

n 

30 

n 

32 

n 

34 

D 

36 

n 

40 

n 

41 

I! 

42 

D 

46 

n 

50 

n 

66 

n 

74 

ii 

75 

ii 

7.J 

n 
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80  n     input  n,:  4  sexadecimal  characters  from  tape 

82  n     punch  n/4  sexadecimal  characters  on  tape 

K5  n     put  Q  +  2~~39into  A 

S5  n     put  Q  into  A 

JO  n     replace  Q  by  logical  product  of  Q  and  N(n) 

F5  n     put  N(n)  +  2~39  into  A 
L0  n     subtract  N(n)  from  A 

LI  n     put  -  N(n)  into  A 

L2  n     subtract  |N(n)  |  from  A 

L3  n     put  -  |N(n)  |  into  A 

L4  n     add  N(n)  to  A 

L5  n     put  N(n)  into  A 

L6  n     add  |N(n)  |  to  A 

L7  n     put  |N(n) |  into  A 

A  number  of  92  output  orders  for  printing  special  teletype 
characters  are  available,  a  few  will  be  listed    All  letters  and 
punctuation  marks   can  be  printed 

92  963     space 

93  131     carriage  return  and  line  feed 
92  515     delay 

92  259  letters  shift 

92  707  numbers  shift 

92  387  A  or  ) 
etc. 


APPENDIX    III 
SERIES    INVOLVING   SPHERICAL   BESSEL    FUNCTIONS 

The   two   series    to   be   calculated   are 

oo  *■ 

00 

K{^^\  (in-2) 

where  the  a  ,  b  are  constants  that  were  obtained  in  Appendix  I,  and  are 
different  for  each  n  and  m.   Because  the  second  calculation  is  almost 
identical  to  the  first,  only  the  latter  will  be  described.   If  the 
spherical  Bessel  function  of  the  first  kind  is  denoted  by  j   and  that  of 
the  second  kind  by  n  ,  then  the  Hankel  is 

hn(2)(kr)  =  jn(kr)  -  j  nn(kr)  (III-3) 

where  the  j  without  subscript  is  the  operator  of  the  Argand  diagram. 
The  real  and  imaginary  parts  of  the  spherical  Bessel  function  fraction 
are 


Re 


nh 


na 


ir,(ka    i_    i    -  n    i_)    +   n.(ka   n„    i    -  n   n_ ) 

J-Hl JlLrl Jjil nl D-L_ ni  (HI-4) 


[rhn]M  (ka  jn_1   -  n   jn)2    +    (ka   nn_1    -  n   nn) 


lm    f    "hna        . ka<Jn   "n-1   ~  Jn-1   "n> ff  (in..5) 

\  LrhnJa  /        ^ka  Jn-1    ~  n   Jn^      +    (ka   nn-l   "  n   nn' 

All    spherical    Bessel     functions   are   evaluated   at    the   argument   ka. 
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Scaling  is  a  serious  problem  here  as  the  spherical  Bessel  functions 
cover  a  wide  range  of  values.   To  obviate  this,  use  has  been  made  of  a 
floating  decimal  library  program.   This  routine  is  an  interpretive  code, 
i  e  ,  the  routine  in  effect  forms  a  new  machine  that  handles  numbers  in 
floating  decimal  form.   It  has  a  floating  accumulator  (in  the  memory)  and 
a  new  order  code   The  problem  is  coded  as  usual  except  that  the  new 
order  code  is  used  and  control  is  transferred  to  the  subroutine   The 
subroutine  interprets  each  new  order  and  follows  the  proper  sequence  of 
regular  orders  to  convert  from  floating  decimal  to  binary,  add,  subtract, 
scale,  reconvert,  etc.   The  limitation  of  these  codes  is  their  intrinsic 
slowness,  by  a  factor  of  about  20  (50  for  addition  and  10  for  division), 
Coding,  on  the  other  hand,  is  greatly  simplified  because  the  interpretive 
subroutine  has  internal  provisions  for  handling  red  tape  or  bookkeeping 
processes. 

No  library  routines  to  calculate  spherical  Bessel  functions  were 
available   thus,  values  were  read  from  tables  and  punched  on  tape  to  be 
read  into  the  machine   The  first  step  of  the  calculation,  after  the 
numbers  have  been  read  in  and  stored,  is  to  compute  Eqs   III-4  and 
III  -  5   These  are  printed  as  they  are  calculated   After  this  has  occurred 
for  n  up  to  ten,  the  terms  of  the  series  are  formed  and  the  series  summed 
for  each  m   These  results  are  also  printed  out   A  block  diagram  is 
shown  in  Fig  III- 1  and  a  typical  output  in  Table  III - 1 
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+999341 

+00 

+163271 

+00 

+995193 

+00 

+344473 

+00 

+980108 

+00 

+565257 

+00 

+927879 

+00 

+866987 

+00 

+724807 

+00 

+132496 

+01 

-893408 

-01 

+181255 

+01 

-114535 

+01 

+985087 

+00 

-817478 

+00 

+175458 

+00 

-486961 

+00 

+228330 

-01 

-325720 

+00 

+249976 

-02 

-270211 

+00 

+260340 

-03 

+00 

-64 

+895596 

+00 

+506870 

+00 

+20 

+01 

+109263 

+00 

+349327 

+00 

+40 

+01 

-753936 

-01 

+196348 

+00 

+60 

+01 

-109853 

+00 

+  560569 

-01 

+80 

+01 

-357072 

-01 

+892486 

-03 

+10 

+02 

-126242 

-01 

+751196 

-05 

+10 

+01 

+354415 

-00 

+448749 

+00 

+30 

+01 

+167808 

-01 

+290973 

+00 

+50 

+01 

-846435 

-01 

+137149 

+00 

+70 

+01 

-599906 

-01 

+816438 

-02 

+90 

+01 

-168156 

-01 

+832884 

-04 

Table 

III- 

.  Typica 

il  Printed 

Results 

APPENDIX  IV 
WATSON'S  TRANSFORMATION 


This  material  represents  an  exploratory  investigation  of  Watson's 
Transformation,  which  is  a  method  for  summing  slowly  convergent  series. 
This  method,  first  used  for  problems  involving  propagation  around  the 
earth,  consists  of  three  separate  steps.   First,  the  series  is  replaced 
by  a  contour  integral  taken  around  a  set  of  simple  poles  on  the  positive 
real  axis.   The  residue  of  the  integrand  at  each  pole  gives  a  term  of 
the  series.   Second,  the  contour  about  the  poles  on  the  real  axis  is 
replaced  by  an  exactly  equivalent  contour  about  a  set  of  complex  order 
poles,  so-called  because  they  occur  for  certain  complex  orders  of  the 
functions  involved  (here  spherical  Bessel  functions),  these  poles  are 
in  the  upper  right  half  plane.   Third,  the  contour  integral  is  replaced 
by  the  sum  of  the  residues  at  the  new  complex  order  poles.   The  resultant 
series  is  much  more  rapidly  convergent. 

The  method  was  used  by  Watson ^4JJ  £n  1917  or  earlier,  and  later  by 
Sommerfeld  ^^^'      However  Bremmer  appears  to  have  been  the  first  to  make 
numerical  use  of  the  technique .  ^o,  44.)   Recently  it  has  been  used  by 
Sensiper  and  Lucke  . '  ^' 

The  series  of  interest  is  S(m,a)  given  by  Eq„  4-14  and  the  simple 
case  of  m  ::  0  may  be  taken    Then 


S(0fi 


n  =  l 


(2n  '  1)  ka  hna  Pn'(0)  Pn'(A) 
2n(n  +  1)  [rhn]a' 


(JV-1) 
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It  will  be  shown  that  this  series  is  exactly  equal  to  the  following  contour 
integral,  taken  around  the  path  of  Fig   IV-1 

S  =  i  /     v   ka  h^  P^(0)  P^(A)  (TV  9) 

J  Ci  2{v  -   X)(v  +  %)   cos  vn    [rhvJ/2]a  UV"*' 

Sec  vtt  has  simple  poles  on  the  real  axis  at  vk   =  rc/2   371/2,  .  .  .  ,  or  v   - 

1/2,  3/2,  5/2,  7/2,...    Next  the  residue  of  the  integrand  at  each  of 

these  must  be  found   The  appropriate  Lemma  from  complex  variable  theory 

is 

Lemma   Let  a  function  G(z)  have  simple  zeroes  at  z0:  zlj    z2/  etc...  and 

let  F(z)  be  finite  at  each  of  these   Then  the  residue  of  F(z)/G(z) 

at  any  of  these  poles  is  simply  F(z0)/G'(z)    , 
The  proof  follows  immediately  from  the  method  for  finding  residues  when 
only  simple  poles  are  involved   Thus 


Residue  of  F(z)/G(z)  =  Lim   [(z  -  z0]  F/G]  =  F(z0)/G'(z) L        (IV-2) 

z-'Zo  '  ° 

at   Zn 


It  can  now  be  observed  that  the  residue  of  Eq  IV-2  times  2rtj  yields  one 
term  of  the  series  (IV-1)  The  sign  change  (-l)v~'2  has  been  incorporated 
into  P^  [cos  (71  -  0)]  ~-   PV^(A) 

Next  the  different  and  unrelated  contour  C2  of  Fig  IV-1  will  be 
considered,  the  integral  around  C2  is  zero   This  path  is  along  the  imagi- 
nary axis  and  around  a  semicircle  of  radius  R,  the  semicircle  has  center  at 
the  origin  and  the  imaginary  axis  is  the  diameter   R  tends  to  infinity 
A  brief  outline  of  the  proof  will  be  given   First   the  integrand  of 
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Figure  U-l        Paths   of    integration 
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Eq   IV-2  must  have  no  singularities  on  C2   This  is  true  because  for  Re 
(argument)  >  0,  both  the  Legendre  function  and  the  spherical  Bessel 
functions  are  entire   Second,  the  integrand  except  for  v   is  even.   From 
the  hypergeometric  function  representation  of  Pv>  Pv_i/  (x)  =  P_j,_y  (x), 
so  the  Legendre  function  here  is  an  even  function  of  v.      Cos  vn   is 
clearly  also  even.   For  the  spherical  Bessel  function  term,  let 

h„  1/  h, 


F(v)  -    ^    =  ^ 

_d_ 

dx 


d-     [xh^]  xhV-3/2  (V    -    Y2)     \lvJ£ 


Now   for   any  v,    ^ly\y2  (x)    =   e'3vK  h^f^   (x)   which  gives 

f(-v) ^ i^V^ H  F(v) 

xh_v_3/2    +    (v   +  K)    h.v.X       e'J^[xhv_3/2   -     (v  -   K)   hv„i/z] 

and  thus  this  factor  is  also  even.   A  recursion  formula  was  used  in  the 
last  step.   All  the  integrand  is  even  except  v   which  makes  it  all  odd 
Contributions  over  paths  A  and  B  will  then  cancel.   From  the  asymptotic 
expansions  for  large  order  of  Legendre  functions  and  spherical  Bessel 
functions,  one  may  easily  show  that  as  the  radius  of  the  semicircle  of 
C2  tends  to  infinity,  the  integrand  tends  to  zero.   Then  by  Jordan's 
Lemma,  the  integral  over  the  semicircle  tends  also  to  zero,   The  integral 
is  now  zero  over  path  C2 

It  was  mentioned  earlier  that  the  spherical  Bessel  function  factor 
in  the  denominator  is  zero  for  certain  complex  orders;  hence  the  inte- 
grand is  singular  at  these  v  Note  that  the  argument  of  the  spherical 
Bessel  function  is  fixed,  and  the  zeroes  (simple)  occur  for  complex 
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order.   Now  the  integral  around  C2  is  zero  so  the  sum  of  residues  on  the 
real  axis  must  be  exactly  the  negative  of  the  sum  of  residues  at  the 
complex  order  poles,   The  transformation  has  now  been  effected,   The 
resulting  series  is 

V"       v  ka  hv  u  PI   i/(0)  P'  i/,(A) 
S  =  -2ti  /   v  n     v  /2      — v  "■      ■.    (IV-4) 

V  2(v  -  Y2)(v   +  %)   cos  vk  $-  {[TU     i/]'} 
vm  3v     v   /2  a 

This  series  is  exact;  no  approximations  have  yet  been  made.   Some  ap- 
proximation will  have  to  be  made  though,  because  at  present  no  means 
exist  for  locating  the  complex  zeroes  and  for  evaluating  the  derivative, 
with  respect  to  v,    at  these  zeroes.   The  single  exception  is  that  of 
large  argument  ka,   If  ka  >>  0,  then  one  of  the  several  asymptotic 
formulas  may  be  used.   Bremmer  has  shown  that  zeroes  occur  only  for 
large  v   which  means  that  an  expansion  for  large  argument  and  order  must 
be  employed   Watson  and  Sommerfeld  used  the  Debye  type  approximation 
for  large  ka  and  v   which  is  obtained  by  the  method  of  Steepest  Descents. 
In  this  method,  the  Hankel  function  is  written  as  an  exponential  contour 
integral  with  the  contour  passing  over  two  saddle  points  (two  must  be 
used  for  v   near  ka).   Next  the  path  is  made  that  of  the  steepest  descent 
so  that  only  the  region  of  the  saddle  top  is  of  importance.,   The  exponent 
is  expanded  into  a  three  term  Taylor's  series  about  the  point,  and  the 
integral  can  readily  be  evaluated   The  range  of  validity  is  |v-ka|<<ka. 


Actually  two  terms  as  the  first  derivative  at  the  point  is  zero;  this 
i   in  fact  the  condition  for  the  location  of  the  saddle  point. 
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Use  of  this  approach  demonstrated  that  the  complex  order  zeroes  of 
importance  occur  for  v  very  near  ka  For  v  -  ka,  the  second  derivative 
in  the  Taylor  expansion  is  zero,  thus  for  v  near  ka  the  Debye  result  is 
poor.  This  dilemma  was  resolved  by  Bremmer  who  recognized  that  a  four 
term  Taylor  series  should  be  used,  leading  to  a  third  order  approxi- 
mation The  resulting  contour  integral  cannot  however  be  simply  evalu- 
ated  It  is  an  integral  of  the  Airy  type,  a  simple  example  of  which  is 


/. 


I°°cos  (e3  -  xe)  de  =  ^K!/   '2^ 

*  3       /j 


3nT3~ 


(IV- 5) 


Use  of  a  similar  integral  reduces  the  form  to  a  Hankel  function  of  order 
1/3.  This  approximation  is  called  the  Hankel  approximation.  Applied  to 
the  spherical  Bessel  function  factor, 


[rlv*]«      i    ,  ■ 


kah7.  i/ 
x   =   ka. 


x      x 

-2"        —•. 
V  V 


1 


lA 


X        X 

—  I    — —T 

v      v' 


1 


% 


■j(2n/3)   Jl 


H^>    (3) 


ntfj  O) 


(IV-6) 


Bremmer  has  shown  that  near  the  zeroes,  the  second  term  is  slowly  varying 
and  the  location  of  them  is  thereby  controlled  by  the  factor  Hkr'  ((3). 

'73 

The  first  term  of  course  is  very  small    In  the  notation  above, 


v  -  x 


y. 


3  =   1/3(2t)3/^  e-J<3*/2) 


(IV-7) 


(IV-8) 


.(2) 


(46) 


The  H2/       factor   is    zero    for   the    following  values,  v*o;    where 

Hi)  m 


(IV-9) 
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m 


is  0  am  is  0.6855  xm  is  0.808  eJ(n/3) 

1  3.903  2.577  eJ(7l/3) 

2  7.055  3.824  eJ(Tl/3) 

3  10.200  4.892  eJ(n/3) 


4 


13.345  5.851  eJ(n/3) 


5  16.488  6.737  eJ(rc/3) 

Now  that  the  zeroes  have  been  located,  the  residue  at  each  must  be 
evaluated.   Recourse  must  here  be  made  to  the  Debye  approximation;  now 
that  the  more  accurate  Hankel  approximation  has  been  used  to  locate 
the  zeroes,  the  derivative  may  be  evaluated  by  the  less  accurate 
formulas.   The  result  is 

9.    j     hK-%K     \  a    -77-^"  (IV- 10) 

dv    |ka  h„^(ka)J  <ka>/s    '  UV  i0) 

The   spherical    harmonic   P  can   be   dealt  with    from   the   simple   asymptotic 
form: 

d<        r          <          mi        ^  cos    [y(rc  -  9)   -  n/4]  ,Tlfll* 

?l_y7    [cos    (rt  -   0)]    ~  .  (IV-11) 

•FH  sin   9 

For   positive    Im    (v)  t    in    the   cosine    term   the   second  exponential    is   pre- 
dominant.     Similarly    for   cos    vk. 

P'   „   (A)  -   rZ  e~J(vn/2)    -  j(n/4)    +  jvA 

cos   un  ~  %e'iv1t   .  (IV- 12) 
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So 

COS  VTC  K 

The  negative  real  part  of  the  exponent  means  that  the  values  of  v   near  the 
origin  are  the  most  important 

These  are  inserted  into  Eq   IV-4  to  give 


^~~       P'  i/(0)  P'  i/(A)          n  \2/3  ^~~   eJ2>rn 
S  =  -Tt  )   *V^u;  vv_W*) ^  _.  (ka)    )       f .     (iv-14) 


for  the  final  series   Numerical  values  are  easily  obtained,  using  six 
terms  and  ka  =  6  28  gives 

S  =  -3  2  -  j  16 
which  is  not  close  to  the  result  obtained  by  summing  the  original  series, 
S  -  -4  6  +  j  32   Due  to  the  small  gap  angle,  the  terms  here  do  not  de- 
crease rapidly,  the  sixth  term  is  about  five  percent  of  the  total   This 
is  a  considerable  improvement  over  the  original  series  though,  as  hundreds 
of  terms  were  needed  there   The  reason  for  poor  agreement  is  the  value  of 
ka  of  6  28   This  was  the  largest  diameter  sphere  of  interest  here   but 
the  asymptotic  formulas  are  just  not  accurate  for  this  small  an  argument, 
especially  when  the  zero  point  is  of  importance   This  points  up  the 
limitations  and  usefulness  of  the  method 

The  physical  interpretation  of  Watson's  Transformation  is  of  inter- 
est  Papas ^4''  has  shown  that  rigorous  solutions  to  the  Maxwell  equations 
and  the  wave  equation  can  be  constructed  from  complex  order  wave  functions 
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of  the  type 

h!,2)(kr)  Pv   (cos  9) 

m        m 

where  the  vm   are  the  complex  order  roots  of 

f  [rhv  (kr)]a  -  0 

or     m 

These  solutions  satisfy  the  Sommerfeld  radiation  condition  as  they  must. 

(2) 
The  complex  wave  functions  are  orthogonal  in  an  unusual  way  -  the  hv      (ka) 

m 

are  orthogonal  over  the  various  vm    !   The  Pv(cos  9)  are  no  longer  orthogo- 
nal.  These  complex  order  wave  functions  are  just  the  terms  of  the  residue 
series  of  the  Watson  Transformation. 
Whereas  the  classical  solution 

hn2)(kr)  Pn(cos  9)    n  =  1,  2,  3,  .  .  . 

represents  a  travelling  wave  in  the  r  direction  and  a  standing  wave  in 
the  9  direction,  the  complex  order  solution  represents  a  travelling  wave 
in  the  9  direction  (with  Green's  function  source  at  9  -  0),  and  a  wave 
going  in  the  r  direction  that  is  reflected  at  infinity;  in  other  words  a 
standing  wave  in  the  r  direction  at  infinity.   In  some  physical  problems, 

the  reflection  coefficient  at  infinity  may  be  zero,   Note  that  only 

(2) 
ttjj   (kr)  can  be  used  here  in  the  classical  solution  as  jn(kr)  and  nn(kr) 

do  not  satisfy  the  radiation  condition. 

To  avoid  the  circuitous  path  of  Watson,  a  problem  might  be  set  up 

directly  ID  complex  order  wave  functions    A  worthy  project  might  be  the 

calculation  of  the  zeroes  u      for  various  ka 


APPENDIX    V 


THE    Zm    CALCULATION 


From  Eq-    4-13    it   is    seen    that    Zm    is    a    simple    function    of 
S       S    ,    Jm,    Jm'        The   program    to    perform    the    calculation    uses    the 
floating  decimal   subroutine,     first    reading    in    the    S    functions, 
then   the  J   functions        After    the    results    are    printed    out,    a    new 
set   of  functions   may  be    read    in    and    the    calculation    repeated. 
Several   tables   of  Bessel    functions    were    found    usef  ul  ^°  '  ^'  »  **"  '  ^  ' 
Results   are   tabulated   in    Tables    V-l    and    V- 2 


0 


0 

16 

-2184-02 

-3718+01 

-3319+02 

-6711+03 

0 

31 

-5677+01 

-1910+01 

-8851   01 

-9250+02 

0 

47 

-2380-01 

-1240+01 

-3927-01 

-2684-02 

0 

63 

-1277-01 

-9025-00 

-2155+01 

-1112+02 

0 

78 

-8193+00 

-7042-00 

-1412+01 

-5900+01 

1 

05 

-4842-00 

-4798-00 

-7868-00 

-2438-01 

1 

25 

-4003+00 

-3733+00 

-5977-00 

-1481+01 

1 

s-7 

-3586+00 

-2661-00 

-3583+00 

-7720+00 

1 

88 

-1340-01 

-2165+00 

-2507-00 

-4677-00 

2 

05 

+7602-00 

-2063+00 

-2110+00 

-3694+00 

2 

20 

+8499+00 

-2045-00 

-1830+00 

-3056+00 

2 

40 

+5164  00 

-2100+00 

-1558+00 

-2426-00 

2 

52 

+3822+00 

-2172+00 

-1432+00 

-2135^00 

2 

83 

+2068+00 

-2529+00 

-1225+00 

-158400 

3 

14 

+1384+00 

-3480J00 

-1151+00 

-1225+00 

3 

50 

+1053+00 

-5448+00 

-1156+00 

-9666-01 

3 

76 

+9411-01 

+5501-00 

-1216-00 

-8592-01 

4 

71 

+9973-01 

+9451-01 

-3869-00 

-7845-01 

5 

34 

-2790+00 

+6844-01 

+  1909^00 

-7177-01 

5 

80 

-1359+00 

+6883-01 

+9203-01 

-2081+00 

6 

28 

-6249-01 

+9459-01 

+6265-01 

-6273+00 

Table 

V-l        Real 

Part  of  Zm 
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m  4           b           b  I 

ka 

0,16  -2164+05  -9644+06 

0.31  -1539+04  -3540+05  -1043+07 

0.47  -2937+03  -4449+04  -8654+05  -2427+06 

0.63  -9112+02  -1031+04  -1493+05  -2639+06 

0.78  -3906+02  -3581+03  -4170+04  -5947+05 

1.05  -1209+02  -8197+02  -7109+03  -7528+04 

1.25  -6114+01  -3477+02  -2532+03  -2249+04 

1  57  -2534+01  -1145-02  -6621+02  -4674+03 

1,88  -1266+01  -4817+01  -2318+02  1363+03 

2.05  -9280+00  -3191+01  -1407+02  -7570+02 

2,20  -7076+00  -2205+01  -9384+01  -4699+02 

2.40  -5236+00  -1529+01  -5574+01  -2621+02 

2.52  -4405+00  -1223+01  -4346+01  -1894+02 

2.83  -2949+00  -7243+00  -2281+01  -8800+01 

3.14  -2084+00  -4604+00  -1295+00  -4476+01 

3.50  -1488-00  -2239+00  -4914+00  -1374+00 

3.76  -1179+00  -2181+00  -5039+00  -1433+01 

4.71  -6640+01  -9387-01  -1705+00  -3764+00 

5.34  -5828-01  -6317-01  -9908-01  -1901+00 

5  80  -5839+01  -5254-01  -7148-01  -1252+00 

6.28  -6620-01  -4780-01  -5420-01  -8616-01 


ka 


10 


0  16 

- 

- 

- 

0  31 

- 

- 

- 

0  47 

- 

- 

- 

0  63 

- 

- 

- 

0  78 

-1004  07 

- 

- 

1.05 

-9424+0  5 

-1364-07 

- 

1  25 

-2363-0  5 

-2865+06 

- 

1  57 

-3900+04 

-3761-04 

-4113+06 

1  88 

-9378+03 

-7616+04 

-6940+05 

2.05 

-4818+03 

-3546+04 

-2961+05 

2  20 

-2783+03 

-1906*04 

-1510+05 

2  40 

-1419+03 

-8901+03 

-6327+04 

2  52 

-9748-02 

-5814+03 

-3933+04 

2.83 

-4017+02 

-2127+03 

-1277+04 

3.14 

-1833+02 

-8708+02 

-4696+03 

3.50 

-457  5+01 

-1769+02 

-1464+03 

3  76 

-4842  01 

-1902+02 

-8490+02 

4  71 

-1029^01 

-3039*01 

-1063+02 

r>  34 

4301+00 

-1148  01 

-3489+01 

5  BO 

-2567-00 

-6214+00 

-1716*01 

6  28 

-1601 +00 

-3495*00 

-8310+00 

Table  V-l 

Real  Part 

of  Z„  (conto ) 
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0   16 

-5955-01 

+2605-02 

- 

- 

0.31 

1165+00 

+8950-02 

+2728-03 

_ 

0.47 

-1757-00 

+1754-01 

+9565-03 

- 

0  63 

-2341+00 

+2533-01 

+2094-02 

+1244-03 

0,78 

-2872+00 

+3029-01 

+3635-02 

+2841-03 

1,05 

-3811+00 

+2989-01 

+7143-02 

+8605-03 

1.25 

-4586+00 

+2185-01 

+9802-02 

+1587  02 

1   57 

-6665+00 

-1683-02 

+1268-01 

+3192-02 

1,88 

-1171+01 

-2873-01 

+1252-01 

^5043-02 

2.05 

1013+01 

-4112-01 

+1092-01 

+6003-02 

2.20 

-3343+00 

-4955-01 

+8537-02 

+6713-02 

2,40 

-2947-02 

-5767-01 

+4296-02 

+7353-02 

2  52 

+4399-01 

-6185-01 

+1245-02 

+7516-02 

2  83 

+6625-01 

-7773-01 

-6719-02 

+7010-02 

3  14 

+6420-01 

-1359+00 

-1215-01 

+4998-02 

3.50 

+5840-01 

-1173+01 

-1444-01 

+1077-02 

3.76 

+5500-01 

-6163-01 

-1581-01 

1914-02 

4,71 

+8560-01 

+3165-01 

-1788+00 

-5642-02 

5  34 

+5765+00 

+2886-01 

+2520-01 

-6502-02 

5.80 

-5775-01 

+3687-01 

+2080-01 

-5931-01 

6.28 

-3583-01 

+6886-01 

+1791-01 

-1191+00 

m 

4 

5 

6 

7 

ka 

0    16 

_ 

„ 

„ 

=. 

0.31 

- 

- 

- 

- 

0,47 

- 

- 

- 

- 

0,63 

- 

- 

- 

- 

0.78 

- 

- 

- 

1.05 

+7221-04 

- 

- 

- 

1.25 

+1726-03 

- 

- 

- 

1.57 

+5088-03 

+7241-04 

■ 

— 

1.88 

+1042-02 

+1806  03 

- 

- 

2  05 

+1483-02 

+2590-03 

- 

- 

2.20 

+1872-02 

+3490-03 

+6287-04 

- 

2.40 

+2572-02 

+5871-03 

+1183-03 

- 

2.52 

+2986-02 

+7297-03 

+1320-03 

- 

2,83 

+4017-02 

+1229-02 

+2696-03 

+6019-04 

3  14 

+4814-02 

+1860-02 

+4922-03 

+1011-03 

3  50 

+5220-02 

+2167-02 

+5697-03 

+1274-03 

3  76 

+4724-02 

+3218-02 

+1211-02 

+3428-03 

4  71 

-5109-03 

+3520-02 

^2683-02 

+1180-02 

5  34 

-2336-02 

+1403-02 

+305402 

+1971-02 

5  80 

-2443-02 

-4051-03 

+2559-02 

+2437-02 

6.28 

-5356-02 

-1160  02 

+1297-02 

+2589-02 

Table  V-2 

Imaginary 

r  Part  of  1„ 
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10 


2.83 

- 

- 

- 

3.14 

- 

- 

- 

3.50 

+3219-04 

- 

- 

3.76 

+7179-04 

- 

- 

4.71 

+4183-03 

+9496-04 

- 

5.34 

+8261-03 

+2713-03 

+7981-04 

5.80 

+1261-02 

+4916-03 

+1551-03 

6.28 

+1738-02 

+7974-03 

+2545=03 

Table  V-2.  Imaginary  Part  of  Zm  (Cont  ) 


APPENDIX  VI 
THE  ASSOCIATED  LEGENDRE  CODE 

The  ILLIAC  program  to  calculate  associated  Legendre  polynomials  (and 
the  auxiliary  program  to  calculate  the  derivative)  is  arranged  as  a  stand- 
ard closed  subroutine  so  that  it  can  be  used  by  others.   Entry  is  made 
with  the  parameter  n  in  the  link  and  the  angle  0°  *  2~ *■*   (as  a  left  hand 
address)  in  the  accumulator.   The  library  sine-*cosine  and  square  root 
codes  must  be  appended  to  this  code.   Also  1  £  n  ±   10  due  to  scaling 
difficulties.   The  routine  is  vacated  with  em/4  times  the  normalized 
associated  Legendre  polynomials  Pn°(cos  9),  Pni(cos  9),  ....  P  n(cos  9) 
in  the  preset  locations  S3,  1S3,  .  .  .  ,  nS3.   Scaled  this  way,  the  functions 
are  always  less  than  unity. 

To  avoid  singularity  difficulties  at  9   0,  calculation  is  first  made 
of  the  R  functions,  defined  here: 


Pnm(cos  9)  =  sinm  9  Rnm(cos  9)  .  (VI-1) 


Initial  values  for  the  R  function  are  first  calculated,  then  used  in  the 
three  term  recursion  formulas  to  generate  successive  functions.   The 
formulas  are 


r  m  =  y2    |(n  +  lA)(n  -%)...  (Yi)  (VI-2) 


>n-l 


Hg-1  *M   Bn„!  cos 


(n  +  %)(n   -  X)...(X)n       (yi_3) 
2n!  U1  6) 
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R-2   .    £         j     (m-_l)    pos   9Bg-l        _  l/2  sin2    0     r    (n+m)(n-m+U_  or 
^  m-2  M(n  +m-l)(n-m  +  2)  V  (n  +  m  -  1)  (n  -  m  +  2)     ta 

(VI-4) 

After  the  set  of  R' s  has  been  stored,  the  P's  are  found  by  multiplying 

by  sinm  9.   These  new  values  then  replace  the  R' s  in  the  store. 

The  auxilliary  program  calculates  the  derivative  functions,  times 

em- 1^2  •   Utilizing  the  Pnm  code,  it  permits  formation  of  the  set  of  R's, 

then  finds  the  derivative  from  well  known  recursion  formulas  involving 

the  Pnm,  finally  returns  to  the  original  subroutine  and  stores  the  Pnm. 

The  equations  involved  are 


gp  o  

-g-  =   -2~5  sin    6  ^n(n   +   1)    R^ 


39 


^  =   2-4  sxn-1 
39 


-^ —  nI  (n    +  m)  (n 
em-l 


A  block   diagram  will    be    found   in    the  next   pages   along  with   the    codes 
themselves . 

Normalized  Associated  Legendre  Polynomial   Code 


0  40   F 
S5  IF 

1  L4  L 

42   98L     link 

2  10   20F 
42   102L 

3  22   3L 
51   F 

4  00   12F 
75   1061. 

5  00   IF     V2  9 
50  5L 


6  26  112L     to  sin-cos 
40  107L 

7  L5  F 
^40  108L 

8  L5  105L 
L4  103L 

9  42  91L     set  add. 
42  93L 

10  L5  105L 
L4  102L 

11  L0  103L 
42  35L 
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12 

42 

65L 

L0 

103L 

13 

42 

46L 

42 

73L 

14 

LO 

103L 

42 

79L 

15 

L5 

102L 

42 

30L 

16 

L5 

103L 

40 

IF 

17^ 

40 

2F 

"51 

2F 

18 

75 

IF 

00 

39F 

19 

40 

2F 

L5 

IF 

20 

L4 

40 

10  4L 
IF 

21 

10 

IF 

LO 

102L 

22 

36 

23L 
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IF 
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31 

40 

IF 

L5 

2F 

32 

66 

IF 

41 

F 

33 

S5 

F 

50 

33L 

34 

26 

142L 

40 

L11L 

35 

10 

•  I 

40 

(nS3 

His   +   3/ 


21 


38 


r(K) 


n!2 


39 


DF 


n!2 


n-36 


to   square    root 


%\\r 


36  L5  102L 

L0  103L 

37   22  37L 

50  37L 

38  26  142L 

%  U/2 

00  13F 

39  40  IF 

51  IF 

40  75  107L 

00  IF 

41   40   IF 

51   IF 

42  7J  111L 

40   IF 

43  L5  104L 

L0   102L 

44  36  46L 

En~l 

L5  IF 

45  00  IF 

40   IF 

4<Tl5  IF 

40   OF 

^  ! 

47  L5  104L 

L0   102L 

48  36  88L 

is  n   -    1? 

L5  102L 

49  L0   104L 

set   i 

40  111L 

~l5"L5  102L 

L0  111L 

51  40   109L 

(n-m  +  2) 

L0   103L 

52  40   110L 

(n-m+  1) 

L5  102L 

53  L4  111L 

40  IF 

(n  +  m) 

54  LO  103L 

40  2F 

( n   +  m  '  — 

55  51  2F 

75  109L 

56  00  39F 

40   100L 

57   51   110L 

75  IF 

58  00  37 F 

66  109L 

59  S5  F 

50   59L 

1) 


143 


60  26  142L 
40  2F 

61  51   108L 
75   108L 

62  00   IF 
40   IF 

63  51   IF 
75  2F 

64  00  IF 
40  2F 

65  50  2F 

7J   (mS3)F 

66  40   110L 
51   111L 

67  75  111L 
00  37F 

68  66  109L 
S5  F 

69  22  69L 
50  69L 

70  26  142L 
40  F 

71  51   F 

75  107L 

72  00   IF 
40  2F 

73  50  2F 
7J   OF 

74  00  IF 
L0   110L 

75  40   IF 
L5  103L 

76  LO   111L 
32  78L 

77  L5  IF 
00   IF 

78  40   IF 
~~^5  IF 

79  22  79L 
40   OF 

80  L5  103L 
LO   111L 

81  36  88L 
L5  111L 

82  LO  103L 
40  111L 

83  L5  65L 
LO  103L 


;m  -2 


KRK 


-2 


is  1 


1? 


dec   counter 


84  42  65L   dec.  add. 
L5  73L 

85  LO  103L 
42  73L 

86  L5  79L 
LO  10  3L 

87  42  79L 
_^26  50L 

88  L5  103L 

40  F     set  counter 

89  L5  108L 
40  IF 

90  L5  S3 
00  6F 

91  40  S3 

51  (1S3)F 

92  7J  IF 
00  7F 

93  22  93L 

40  (1S3)F 

94  51  IF 
7J  108L 

95  00  IF 

40  IF  lA   sinm+10 

96  L5  F 
L4  103L 

97  40  F    advance  counter 
LO  10 2L 

98  22  98L 

32  OF  exit 

99  L5  91L 
L4  103L 

100  42  91L 
42  93L 

101  22  91L 
OF  F 

102  00  F 

00  (n  +  1)F 

103  00  F 
00  IF 

104  00  F 
00  2F 

105  00  S3 
00  S3 

106  40  F00      )  8_rt 
058505360000J7  45 

107  00  F 
00  F 
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108 

00 

F 

00 

F 

109 

00 

F 

00 

F 

110 

00 

F 

00 

F 

111 

00 

F 

00 

F 

code   66  at 

112L 

code   2'. 

.   at 

142L 

Derivative  Code 


0   40 

F 

S5 

IF 

1  L4 

L 

42 

153L 

2  46 

51L 

L0 

53L 

3  46 

52L 

L5 

23L 

4  42 

56L 

L5 

7L 

5  46 

103L 

alter  NALP  c 

46 

136L 

6  L5 

52L 

46 

8L 

plant  n 

7   50 

10L 

L5 

F 

8   50 

(rOF 

50 

8L 

9  26 

55L 

enter  NALP 

_26 

10L 

10  L5 

160L 

46 

33L 

set   add. 

11   L4 

53L 

46 

24L 

12  L5 

54L 

L4 

53L 

13  46 

41L 

49 

166L 

set  XA  sinm"l 

H  51 

51L 

75 

52L 

15  00 

30  F 

n(n    s   1)2"8 

50 

L5L 

](,  26 

197L 

40 

IF 

17 

50 

IF 

7J 

163L 

18 

40 

IF 

51 

IF 

19 

71 

1S3 

00 

6F 

20 

40 

S4       P] 

L5 

53L 

21 

40 

165L     ! 

to50 

165L 

22 

00 

15F 

7J 

162L 

23 

40 

IF         . 

50 

IF 

24 

75 

(1S3)F 

00 

IF 

25 

40 

164L 

L5 

165L 

26 

L4 

52  L 

40 

IF 

27 

L5 

51L 

L0 

165L 

28 

40 

2F         1 

51 

2F 

29 

75 

IF 

00 

30F 

30 

22 

30L 

50 

30L 

31 

26 

197L 

40 

IF 

32 

50 

IF 

26 

33L 

33 

7J 

(S3)F 

40 

IF 

34 

L5 

165L 

L0 

53L 

35 

L0 

53L)     , 
37LJ 

32 

36 

L5 

IF 

00 

IF 

37^ 

40 

IF 

L5 

IF 

38 

L0 

164L 

40 

IF 

39 

51 

IF 

75 

166L 

40 

00 

7F 

26 

41L 

set  m    =    1 


cos 


(n    +  m) 


(n   -  m   +   1) 


145 


41  40  (1S4)F    Pnm 
L5  165L 

42  L4  53L 

40  165L      inc.  m 

43  L0  51L 

36  143L      reenter  NALP 

44  51  166L 

75  163L     Y2   sinm  9 

45  00  IF 
40  166L 

46  L5  33L 
L4  53L 

47  46  33L       advance  add. 
L4  53L 

48  46  24L 
L5  41L 

49  L4  53L 
46  41 L 

50  22  21L 
00  F 

51  00  (n  +  1)F 
00  F 


52  00 

(n)F 

00 

F 

53  00 

IF 

00 

F 

54  00 

S4 

00 

F 

NALP  i 

it  (p 

+ 

55) 

this  < 

:ode  at 

(P) 

code  ( 

56  at 

(P 

+ 

167) 

code  21  at 

(P 

+ 

197) 

Both  function  and  derivative  are  normalized  by 


I    2(n  +  m)! 
J(2n  +  l)(n  -  m)! 


ERRATA 


Page  45.   Eq.  5-4  should  read 


jni  H02)(kri) 

2rtcH0o 


with   Hqo  = 


3H0(2)(kri) 
Bkri 


Page  68.   The  first  sentence  after  Eq  6-2  should  read: 

''For  small  c  (with  respect  to  \)  the  behavior  of  Z0o  will  be 
inspected. " 


Page  69.   The  fifth  sentence  should  read: 

"Using  Eq.  6-3  to  get  ZCo  in  terms  of  Bessel  functions  and 
S's 


,     s-UtNp  -  Mi) 


'0  0 


S'Jq  +  J5 


(6  9) 


Page  69.   The  second  equation  from  the  bottom  of  the  page  should  read 

_  ,   Sq  (4  +  ^k2c2Na) 


'00 


rckc(2S0-  1) 


and  the  last  equation  should  read 


Zoo  ~  4Sc/7ikc(4Sc  -  1) 


Page  70.   The  second  part  of  Eq  6-10  should  read 


ZoG   : 


Tike 


